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THE USE OF HIGHER-ORDER DIFFERENCE METHODS IN BEAM 
VIBRATION ANALYSIS* 


By Donald T. Greenwood** 


SUMMARY 


Simple and higher -order difference methods for the solution for the natural 
frequencies of vibration of a uniform beam are compared. The same basic higher- 
order method is used throughout for the interior cells , but three different 
methods of boundary -condition representation are given. 

Tables and graphs of the error in mode frequencies , as compared with a 
continuous beam, are given for the various methods as a function of the number 
of cells. It is concluded that higher-order methods improve accuracy for a 
given number of cells, with essentially no change in the quantity of com- 
puting equipment required. 
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ministered through The University of Michigan's Office of Research Adminstra- 
tion, Ann Arbor . 
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I . INTRODUCTION 


One of the most commonly encountered problem; of structural dynamics is 
that of the lateral vibrations of beams . With th 3 advent on a large scale of 
analog and digital computers of greatly improved speed and computational ef- 
ficiency, it is worthwhile to take another look a: the methods available for 
performing these computations. This report is principally concerned with im- 
provements in analog computer methods of analysis of the be am -vibration prob- 
lem. In particular, it is concerned with methods applicable to the electronic 
differential analyzer. Nevertheless, it will be seen that the matrix formula- 
tion of the problem and most of the results are directly applicable to digital 
computation as w r ell . 

Until quite recently, all attempts to use amlog computers in the solu- 
tion of partial differential equations by difference techniques were re- 
stricted to simple methods with accuracies of second order in the cell size 
at best. In 195b, Fisher^ suggested the use of higher-order differences in 
the solution of these problems on the analog computer and indicated some of 
the advantages to be expected. However, he was not concerned in his examples 
with the beam equation and did not study the problem of boundary-condition 
representation using higher-order methods . 

This report considers in detail the adequate representation of boundary 
conditions when higher-order methods are used, and presents the results of 
rather extensive computations which enable one to obtain a feel for advantages 
and disadvantages of each method. The analysis i 1 the main body of the re- 
port is concerned with uniform beams. The mode -frequency errors for a tapered 
cantilever beam are considered in Appendix B and comparisons are made with 
the continuous beam. 

No analysis of mode -shape errors has been included. However, spot checks 
have shown these errors to be comparable in magni dude with the frequency errors . 

The author wishes to acknowledge the contributions and helpful sugges- 
tions of R. M. Howe of the Department of Aeronautical and Astronautical Engi- 
neering, particularly with respect to the section on errors. Thanks are also 
due C. K. Shah, who did most of the numeric al computation, and J. W. Thatcher, 
w r ho aided in the digital computer programming. 
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II. SYMBOLS 

mode -amplitude matrix 

constant coefficients 

-[T]- 1 [K] 

difference matrices 

an element of matrix B 

constant coefficients 

capacitance, constraint matrix 

constants 

difference matrix 

Young’s modulus 

constraint force 

lateral force per unit length 

constant 

c o ns t ant 

cell size 

moment of inertia of cross section, unit matrix 
mode number 

system matrix [mj [k ] 
stiffness matrix 
bending stiffness of a cell 
inductance* total length of beam 



bending moment 
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M 

m mass, mass matrix 

N total number of cells 

n station designation, number of degrees of freedom 
R resistance 

5 S = 0 is the characteristic equation 

T - [K] [B ] 1 . Also the time -dependent factor in the beam deflection 
t time 

V shear force 

x position along the beam 

Y deflection amplitude 

y lateral deflection 

P eigenvalue for a continuous beam of unit .^ength 
A a small perturbation of the designated quantity 
e frequency error 

£ dimensionless lateral deflection 

T| the square of a mode frequency 

9 beam slope 

A eigenvalue for beam of length N 

£ dimensionless position along beam 

p mass per unit length 

T dimensionless time 


p owo -=* 



eigenvalue for difference beam of unit length 
indices referring to row and column, respectively 
reference value 

transpose of the designated matrix 
computed in the modal (y) coordinate system 
approximately equal to 

amplitude of a sinusoidally varying quantity , determinant of a matrix 
row matrix, or usual parenthesis notation 
rectangular matrix, or usual brackets notation 
column matrix, or usual braces notation 
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III. CONTINUOUS BEJ.MS 


A. THE BEAM EQUATION 

We will restrict ourselves to the study of the lateral vibrations of a 
thin beam where linearity is assumed. The partial differential equation des- 
cribing the motion is 


where 


a 2 

ax 2 



+ 



f(x,t) 


x = position along team 
y = lateral deflection 

t = time 

El = bending stiffness 

p = mass per unit length 
f = lateral force per unit length 


( 1 ) 


It can be seen that Eq. (l) is of fourth orler in x and second order in 
t # For our purposes it is often convenient to write four equations of first 
order in x. 
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= & 
dx 

(2) 

M 

= EI ^ 
ox 

(3) 

V 

dM 

W 

dx 

p 

P St 2 

- - ~ * r(*,i.) 

OX 

(5) 


where 

0 = slope of beam 

M = bending mome nts 
V = shear force 


•f^o>vo ti 
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This set of partial differential equations is equivalent to Eq. (l). We will 
see that the analog computer approach to the problem is based upon this latter 
set of equations rather than the original fourth-order equation. 


B. SOLUTIONS OF THE BEAM EQUATION 

1. General Solution . — In general, the parameters p and El are functions 
of x. Initially, however, we will assume that each is constant. Furthermore, 
we will assume that units are chosen such that the coefficients are unity. 

Our principal interest is in the solution of the homogeneous equation, so let 
us set f(x,t) =0. Then Eq. (l) becomes 


+ fbz = o (6) 

5x 4 bt 2 


It may be seen that essentially the same equation can be obtained by trans- 
forming Eq. (l) to dimensionless form and again setting f(x,t) = 0. Let 




y 

L 


£ 


X 

L 



where L is the length of the beam. Then 


an 

dr 1 + dr 2 


(7) 


Equation (6) can be solved by the method of separation of variables. 
The solution is of the form 


2,3 


y 


Y(x) T(t) 


( 8 ) 
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The Free -Free Beam. — For 


the case of a free-free beam. 


conditions are 


the boundary 


M(0) 

= M(l) = 0 

(24) 

v(o) 

= M(l) = 0 

(25) 


Substituting into Eq. (9) > we obtain 


-Ci + C 3 = 0 (26) 

-c 2 + c 4 = o ( 27 ) 


^osh %/p - cos i + (sinh. \/p - sin n/p)c 2 = 0 
|sinh s/p + sin %/^)ci + ^cosh %/p - cos %/p)c 2 


giving, as before, 


22 

Ci 


cos >/"ft - cosh ‘/g 
sinh \Tp - sin T$ 


( 28 ) 

(29) 


(30) 


The eigenvalue equation is the same as for the built-in case, namely, 


1 


cos sffi cosh 


0 


(31) 
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The value p = 0 has physical significance in this case hut since the mode in- 
volves no elastic deformation, it is omitted. The remaining eigenvalues are 
identical with those for the built-in case. 


The amplitude of the motion is 


5 9 


YU) • 20.(1 ...) ,2C 2 ( Je* 

4 .' 8: 5* 9- 

giving the approximate amplitudes 

Y(x) S 2 Cl (i - Vp x + 1*. . ht + .N 

\ 4: 5: 8: 9: / 

I© 1 = 2 + + + 

\ 3: 41 71 8: J 

jip/fef, - Ad + . Ad + .4\ 

\2i 3: 6: 7: y 

‘ - !r ♦ £ 


|M I = 2Cj 

|v | = 2Ca 


(32) 


(33) 

(3M 

(35) 

(36) 


If we compare the solution for the free -free beam with that for a beam 

built-in at both ends, we note some important similarities. We have seen that 

the nonzero eigenvalues are identical. Beyond that, the roles of the deflec- 
tion y and the bending moment M are interchanged in the two cases as may be 

seen by comparing Eqs. (20) and (35) and also Eqs. (22) and (33)* Similarly, 

the roles of the slope 9 and the shear force V are interchanged. This sim- 
ilarity will hold even for the finite difference solution which will be ob- 
tained later. Therefore it will not be necessary to carry out separate calcu- 
lations for built-in and free -free beams. 

4. The Cantilever Beam. — The boundary conditions in this case are 


y(o) = 0 (37) 

9(0) = 0 (38) 

M(i) = 0 (39) 

v(i) = 0 (40) 


Applying these boundary conditions to the general solution given by Eq. (9)> 
we obtain 


Ci + C 3 = 0 


c 2 + c 4 


0 


(41) 

(42) 
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^os s [ft + cosh sTfi^Cx + ^si n ^P + siiii C 2 = 0 
(sinh s/p - sin 'F&Jcx + |cos + cof;h N^pjca = 0 
From Eqs. (43) and. (44', we obtain the eigenvalue equation 


(43) 

(44) 


1 + cos 


\f p cosh n/p 


0 


(45) 


and also 


C2 

Ci 


cos Vp + cosh n/p 
sin "J p + sinh “/p 


(46) 


the ratio being approximately -1 for all modes higher than the first. 
The first five eigenvalues are 


Pi = 3.516015 

= 22 .03449 

P 3 = 61.69721 

P 4 = 120.9019 
P 5 = 199.8595 


Note that, except for i = 1, the value of Pj_ in this case in approximately 
equal to Pi_! for the built-in or free-free case. 

The amplitudes are given by the equations 


Y(x) 



( 47 ) 
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:i Vi pi x + & ♦£§£♦...♦ 5^ * ♦ 2 

5: 91 CiV^: 6: 1 


|h| - -ac Ji + Sh. + ♦ ... ..> 

4 : 8 : c x \ 5; 9: y 

|v | = -2Cif32 P 2 * 3 + P 2 * 7 + ... + 9s/i + P x _ + P 4 * 8 + . 

3: 71 c x V 2: 8: y 


Each of these expressions contains the same powers of x as was obtained for 
the corresponding built-in case in Eqs * (20) to (23)- 

5. The Pinned Beam.— The boundary conditions for a pinned beam are 


y(o) 

= y(D = 0 

(51) 

M(0) 

= M(l) = 0 

(52) 

Applying these conditions to the 

Cx = 

solution given in Eq. { 9 ), we obtain 
C 3 = C 4 = 0 

(53) 

and therefore 

Y(x) 

= C 2 sin '/p x 

( 5*0 

where the eigenvalues are 

q ^2^2 

Pi = 1 It 

(i = 1, 2, 3, ...) 

( 55 ) 



The power series form of the solution is 


Y(x) 



p^x 3 + jjx 5 

3 : 51 



( 56 ) 


Also 


|e 


c 2 'JH[i 


fix 2 + 

2 i 4: 


= c 2 p 



x + 


a 

f^x 3 

3 i 



|v 


2 

CsP 2 



+ 


Pxf. 

21 



(57) 

(58) 

(59) 


The eigenvalues in this case are not close tc those calculated in the 
previous cases. In fact, for the higher modes (large i), the eigenvalues tend 
toward positions midway between those for the built-in, free -free, or cant- 
ilever beams . 
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IV. DIFFERENCES OF ORDER h 2 


Previously we have obtained solutions for the partial differential equa- 
tions describing a continuous uniform beam. Unfortunately , the analog com- 
puter cannot solve partial differential equations directly but must solve the 
set of ordinary differential equations which are obtained by finite difference 
approximations to the spatial derivative. Ihus we will obtain solution ampli- 
tudes at a discrete set of points rather than continuously as a function of x. 
At the given points , however, the finite difference solution should closely 
approximate the solution for the continuous beam. 


A. FINITE DIFFERENCE EQUATIONS 

An analog computer approach to the beam-vibration problem ujsing finite ^ 
difference equations of order h 2 has been given by Howe and Howe and others. 

The procedure in each case is based essentially upon Eq. (60). 


6y 

dx 


i 



-y n + yn+i 

h 


hf_ Sfy 
24 6x 3 


n+o 


(60) 


The interval between stations is designated by h. Subscripts refer to the sta- 
tion at which y (or one of its derivatives) is calculated. The last term is 
the first error term and is not included in the computer mechanization. It is 
given to provide an estimate of the accuracy of the approximation. Equation 
(60) can be derived by writing the Taylor expansions about station n+*r for y n 

dyi 


and 


y n +i and solving for 


after eliminating y n+ i. 


Using the approximation of Eq. (6o) and again assuming that El and p are 
unity, the beam equations (2), (3), (4), and (5) can be written in difference 
form: 




Mr 


n+ 2 


-yn + y n +i 



h 

-M 

+ 9 n4 


h 

-M n + M n+1 


(61) 

(62) 


h 


( 65 ) 
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Vn^ 


( 64 ) 


Note that displacements and "bending moments are calculated at integer 
stations whereas the slopes and shear forces are calculated at half-integer 
stations. If one calculates all quantities only at integer stations, the re- 
sult is to double the effective cell length and thereby to quadruple the first 
error term. Therefore we will use "staggered” stations when finite difference 
approximations to first derivatives are used. 

The analog computer circuit representing a typical internal cell is 
shown in Fig. 1. It requires five amplifiers per cell if one reverses signs 
at adjacent cells. 

Another approach to the problem is to approximate the second spatial 
derivative rather than the first derivative. In this case the finite dif- 
ference approximation is 


afy 

ax 2 


n 


yn-i ~ 2 y n + yn+i h 2 <j 4 v 
h 2 12 Bx 4 


n 


( 65 ) 


and the corresponding beam equations solved on tie computer are 


Mn 


y n -i - 2 y n + y n +i 

h 2 


( 66 ) 


■Mn-i + 2M n - M 1+1 


( 67 ) 


The analog computer circuit for this case requires only 3 amplifiers per 
cell and is shown in Fig. 2. This circuit is subject to exactly the same 
finite difference errors as the 5 -amplifier circuit and, in fact, can be ob- 
tained as a direct reduction of that circuit. 


B. PHYSICAL INTERPRETATION 

There are two ways in which one can give physical meaning to finite dif- 
ference approximations to the beam equations. 


owo t* 
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The first method is to consider a finite difference expression as an 
approximation to a spatial derivative in the continuous beam at a given point. 
For example, the right-hand side of Eq. (64) approximates - ^| n > which, in 
accordance with Eq. (5) for a continuous beam, is proportional to the accelera- 
tion y n of a local differential element. 

The second method is to think of the set of finite difference equations 
as representing exactly a lumped physical system whose response approximates 
that of the continuous beam. The calculated values of accelerations, shear 
forces, etc., at the various stations have their real counterparts in the 
lumped physical system. This approach is helpful in gaining physical insight 
into the meaning of the equations, particularly for simple differences. Two 
lumped physical representations of this finite difference approximation to 
the beam equations are shown in Fig. 3* 

In Fig. 3a the beam is approximated by a series of massless levers con- 
nected by pin joints. At each joint a spring produces a moment (corresponding 
to the bending moment) which is proportional to the difference in slope angle 
9 of the adjacent levers. The mass of the beam is lumped into point masses 
at the joints. The deflection y is measured at the joints while the slope 9 
is measured at the midpoint of each lever. 

Figure 3b shows an electrical circuit analog of the lumped mechanical 
system of Fig. 3a. It also is described by finite difference Eqs. (6l) to 
(64). The displacement velocities are represented by voltages and the shear 
forces and bending moments are represented by currents. The lumped bending 
stiffnesses are represented by inductors, while the lumped masses are repre- 
sented by capacitors. The transformers, each of turns ratio 1 to h, perform 
the coordinate transformations relating slopes to deflections, and can be 
considered as the analog of the massless levers. 


C. BOUNDARY CONDITIONS 

We have previously given the equations stating the boundary conditions 
on a continuous beam for the various cases under consideration. We must now 
write corresponding equations for the finite difference beam. First, however, 
we must decide whether to end the beam at a point where lateral deflection 
and bending moment are computed (the "integer” stations), or at a point where 
the slope and shear force are computed (the "half-integer" stations). 

Experience has shown ^ that higher accuracy can be obtained for a given 
amount of computation (or degrees of freedom) if the end occurs at a half- 
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integer station for built-in or free ends. For the pinned case it is best to 
end the beam at an integer station. 

1. Built-In End. —Previously we have seen that the displacement y and 
slope 0 are zero at a built-in end of a beam. Suppose, for example, that the 
beam ends at station n+^. We must approximate the deflection at a half- 
integer station where it is not usually defined. Using the Taylor expansion 
approach and keeping only the first error term, we find that 


hf. cfy 
8 ex 2 


n+i 


( 68 ) 


Therefore y since y n +i and are zero* we see from Eqs . (6l) and (68) that 


= y 


n+i 


(69) 


implying that there is no lateral displacement at the integer station adjacent 
to a built-in end. 

2. Free End. — At a free end we find that tie bending moment M and shear 
force V are zero. Using reasoning similar to that of the previous case, we 
find that for a beam ending at station n+1 


Mn = Mn +1 = 0 


(70) 


3* Pinned End . — For a pinned end at static i n the boundary conditions 
are quite straightforward, namely, 


y n = ° 

(71) 

O 

tl 

(72) 


In case the pinned end occurs at station n+.i, the boundary conditions 


are 


o\vo o 
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y n = 

-y n +i 

(73) 

I 6 

n 

-Mn+i 

(74) 


4. The Use of Symmetry . —For vibration problems of a uniform beam where 
the same boundary condition applies at each end, one can make use of symmetry 
to reduce the required amount of computation. All the natural modes can be 
classified according to whether the deflection curve has even or odd symmetry 
about the midpoint of the beam. Consider, for example, a beam that is built- 
in at both ends. If one arranges the modes in order of increasing frequency, 
then modes 1, 3, 5, etc., will have mode shapes with even symmetry about the 
midpoint, whereas modes 2, 4, 6, etc., will exhibit odd symmetry about the 
midpoint. (See Fig. 4.) In either event, however, one need analyze only one- 
half of the beam since the motion of the other half can be deduced from sym- 
metry. 

For the case of even symmetry of the deflection curve, the bending mo- 
ment curve will also show even symmetry, whereas the slope and shear force 
curves will have odd symmetry. Conversely, for odd symmetry in y and M, there 
will be even symmetry in 9 and V . 

The above symmetry assumptions are exact and therefore the same fre- 
quencies and mode shapes result as for the case where the equations are writ- 
ten for the complete beam. 


D. ANALOG COMPUTER CIRCUITS 

1. Cantilever Beam .— The analog computer circuit for the cantilever beam 
is useful in illustrating built-in and free end conditions. Two circuits are 
shown in Fig. 5- The first circuit requires 3 amplifiers per cell, the sec- 
ond, 5 amplifiers per cell. These are nominal values, however, and amplifiers 
can be saved in representing the end cells because of the boundary conditions. 
Thus it can be seen that the actual numbers of amplifiers required to repre- 
sent a four-cell cantilever beam are 9 and. 15, respectively. Cells can be 
added by adding standard circuits such as those in Fig. 1 or Fig. 2, repre- 
senting interior cells . 

2. Built-In Beam. — The computer circuits for a six-cell beam that is 
built-in at both ends are shown in Fig. 6. Note that only one-half of the 
beam is actually represented on the computer; the motion of the other half is 
inferred from symmetry. Circuits are shown for even or odd symmetry about the 
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center of the span. In each case the center occurs at a half -integer station. 
Circuits for the case where the center occurs at an integer station could 
have heen derived by a similar process. This situation would arise for a 
built-in beam with an odd number of cells, assuming, of course, that the ends 
occur at half-integer stations. 

J. Pinned Beam . — The computer circuits representing a six-cell beam 
with both ends pinned are shown in Pig. 7* The circuit for the case of odd 
symmetry about the center. Fig. 7(b), is particularly simple because of the 
additional symmetry about the quarter points alorg the span. 

The circuits to be used in the higher-order methods that we consider are 
quite similar in general form to those of Figs, ^(b), 6, and 7- However, 
most amplifiers will have four rather than two inputs . 


E. MATRIX FORMULATION OF THE DIFFERENCE EQUATIONS 

In the analysis of specific cases of finite difference beams, it is con- 
venient to write the equations in matrix form. It is particularly important 
to use this formulation if digital computations are to be performed. 

In this report the beam equations are, in general, written as four first- 
order difference equations. Using matrix notation, they are as follows: 


(©} 

= [BxJ 

{y) 

(75) 

CM} 

= [B 2 ] 

(©} 

(76) 

(V} 

= [B 3 ] 

CM} 

(77) 

m 

= [B 4 ] 

(V) 

(78) 

m + 

[K] (y} 

= 0 

(79) 


where 

= -[B 4 ][B 3 ][B 2 ][B: ] 


[K] 


(80) 


-f=~ o\vo o 



U\MJ 


21 , 


Although the B matrices are rectangular in general, the K matrix must be square. 
Hie natural frequencies are obtained from the determinants! equation 


K - A 2 ! | = 0 


(81) 


where X is the natural frequency in rad/sec and I is the unit matrix. For the 
case where the K matrix is n x n, there will be n values of A 2 (roots) cor- 
D responding to the n natural frequencies of the system. 

As an example , consider the case of a six-cell team that is built-in at 
both ends . We will write the equations for modes with even symmetry at the 
center. As before , El and p are assumed to be unity. The matrix equations 
are 


2 

% 

I2J 


“1 C T 


r ■> 

y2 


4 


1 


ya 


V. J 


Mi 

m 3 

v* J 



fv 


al 


2 

< y 
% 



1 

-1 



!y2 




ys 



( 82 ) 


(83) 


( 84 ) 


(85) 


Assuming a unit cell size (h = l) and performing the matrix multiplica- 
tions, we obtain 
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N 

ys 


6 -3“ 


ys 

< 

w + 


1 

. r 

Y 3 


-3 2 


Y 3 



u— 




0 


According to Eq. (8l), the roots are obtained from 


(86) 


(6-A 2 ) -3 

-3 (2 -A 2 ) 


0 


(87) 


giving the characteristic equation 


A 4 - 8a 2 + 3 = 0 


( 88 ) 


The roots are 


A?, 2 = 0.394449, 7.60555 

or 

Ai >2 = 0.628052, 2.75732 


These are the actual frequencies for the moles with even symmetry and 
h = 1. For our purposes , however, it is convenient to assume that the total 
length of the beam is unity. This implies that 1 = l/N, where N is the number 
of cells, and results in multiplying each frequency by N 2 . Denoting this 
normalized frequency by 0 ), we obtain 

CD! = 22.6099 
0)2 = 99-2815 
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The same basic matrix formulation will be used in the computations through- 
out the report. The use of higher -order differences will add to the numerical 
complexity but will not change the form of the equations . 


F . RESULTS 

The results of computations of the normalized frequency parameter cd using 
differences of order h 2 are summarized in Tables II and III. The computations 
are for built-in (or free -free) and cantilever beams for various numbers of 
cells. The finite difference error is plotted against the cell number N in 
Figs. 15 and 14 . These results are essentially the same as those given pre- 
viously by Howe and Howe^ and Michie,^ except that the higher mode fre- 
quencies are also included. 


G- ORTHOGONALITY 

We have seen previously that the equations of motion can be put in the 

form 


(y) + [K] {y} = 0 


(79) 


We can also write 


[m] (y] + [k] (y) = 0 


(89) 


where m and k are the mass and stiffness matrices , respectively. Comparing 
Eqs . (79) and. (89)> we see that 


[K] = [mf 1 [k] 


(90) 


Now, for the case of the first-order difference equations under considera- 
tion here, both the m and k matrices are symmetric. In fact, the m matrix is 
diagonal. It is always possible in these linear bilateral undamped systems 
(whose kinetic and potential energies can be written as positive-definite 



quadratic forms) to find a coordinate transformation which will simultane- 
ously diagonalize the m and k matrices. In this c ise the natural modes of the 
system are said to he orthogonal. 

The orthogonality of the modes is not surprising since we have already 
seen (Fig. 3 ) that these equations represent exactly a lumped mechanical or 
electrical system composed of linear bilateral elements. However, when we 
discuss higher-order difference methods, we will find that the k matrix is 
not always symmetric, in which case the modes are not orthogonal. 
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- V. DIFFERENCES OF ORDER h 4 


We have given a brief review of the calculation of eigenvalues for a uni- 
form beam, using simple finite difference approximations to the spatial deriv- 
atives • In this way we have obtained some physical feeling for the main fea- 
tures of this method and perhaps some insight into the more promising approaches 
to the use of higher-order differences. In the remainder of the report we 
will consider several higher -order methods and indicate some of the advantages 
as well as problems in their use. 


A. CHOICE OF THE BASIC DIFFERENCE EQUATION 

Even before one gets into the questions concerning the representation of 
boundary conditions, there are a number of possible choices one can make con- 
cerning the basic higher-order difference equation to be used. First is the 
question whether the approximation ought to be with respect to the first, sec- 
ond, or possibly the fourth spatial derivative. We can immediately reject 
the latter possibility on the grounds that it is quite sensitive to component 
errors in the analog computer circuit. ^>9 Also there is the disadvantage that 
other variables of interest such as 9, M, and V are not directly available 
even though three amplifiers are required per cell. 

The choice between approximating two second-order equations or four first- 
order equations is more difficult. Let us consider each possibility in turn. 

The most obvious approximation of order h 4 for the second derivative with 
respect to x is given by the central difference equation 


£y 

Bx 2 


~ - y n -2 + l6 y n -i - 30y n + l6y n+1 - y n +g h 4 B 6 y 

12 h 2 90 Bx 6 


(9D 


Using this approximation, the equations set up on the computer are of the form 
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^yn -2 + "^ y n-i " 30y n + l6y n+i - yn+^ 


(92) 


n 12h 2 


^n-2 ' l6M n-! + 3 QM n - l6M Q+1 + 


(93) 
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where unity values are assumed for El and p. The analog computer circuit re- 
quires three amplifiers per cell and up to five i lputs per amplifier. (See 
Fig. 6.) 

6v 

Consider now some of the higher-order approximations to — The usual 
central difference approximation of order h 4 is 


6y 

6x 


n 


1 

12h 



8y + 
J n- 1 


8* 


n+i 



h 4 6 5 y 
30 8x 5 


(94) 


The resulting beam equations are 



where El and p are assumed to be of unit value. Ehe analog computer circuit 
requires nine amplifiers per cell in this case because each of the four output 
variables must be generated with both signs. (Ea:h cell would require 2 
integrators, >3 summers, and 4 inverters.) 

Another possibility is to try a higher-order approximation using stag- 
gered stations, i.e., calculate the derivatives ai the midpoints between sta- 
tions where the function is defined. A central difference approximation of 
this sort is the following: 




2 7y n + 2 7y n +-L - f n + 2 
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This approximation is seen to have a smaller error term than that for 
Eq. (94) . Furthermore, it will be seen that less analog computer equipment 
is required to implement it. 

Using the approximation of Eq. (99 ) > the beam equations are 


0 n+i ^ - 27y n + 2 ^n+i ' ^n+a) 

M n = — (9 o - 276 i + 276 , - 6 i3 ) 

24h v 1 -^ n 4 n+ ? n+ ly 

v »*i ■ ik (9 i ■ 27Mn + 2W "* 1 ■ "»*) 
y » ■ ± * 27v »-i - 27 vi * v »*#) 


where, as before, unit values are assumed for El and p. Note that y and M are 
calculated at integer stations whereas 9 and V are calculated at half-integer 
stations . 

The analog computer circuit representation of Eqs . (100) to (103) re- 
quires 5 amplifiers per cell with up to 4 inputs per amplifier. The circuit 
for a typical cell is shown in Fig. 9* The signs of amplifier outputs alter- 
nate on successive cells. 

In spite of the fact that the circuit of Fig. 9 requires 5 amplifiers per 
cell compared to 3 amplifiers per cell for the circuit of Fig. 8, the computa- 
tions in this report are based on the beam equations as given by Eqs. (100) 
to (103)- The reasons for this choice are (l) the greater accuracy of the 
basic difference equation, (2) greater ease and flexibility in applying bound- 
ary conditions and forcing functions, and ( 3 ) less sensitivity to computer 
component errors . 


B. BOUNDARY CONDITIONS AND RESULTS 

Earlier we saw that the representation of boundary conditions was relatively 
straightforward for the case of difference equations of order h 2 . When one 
uses higher-order differences, the problem of applying boundary conditions be- 
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comes more complex. Many approaches are possible and the accuracy of the re- 
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Thus we find that we can omit y± as a coordinate in our analysis because 
of the constraint expressed in Eq. (105)- This was also true when simple ap- 
proximations of order h 2 were used, but in that case y x was equal to zero. 

In a similar manner, using the approximations of Table I and Eq. ( 99 ) > 
we can write the remaining beam equations for the case of a built-in end. 
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Let us apply Eqs. (106) to (109) for the specific case of a 6 -cell "beam 
that is built-in at both ends. Only the modes with even symmetry about the 
center will be considered. The equations of motion are 



Letting h = 1 and applying Eqs. (80) and (8l), we obtain 
I (5.22099-A 2 ) -2.77553 I 


( 110 ) 


( 111 ) 


( 112 ) 


(113) 


-4.29387 


(2. 85815 -A 2 ) 


0 


(114) 
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giving the characteristic equation 


X 4 - 8.07914X 2 + 5- 00^-62 = 0 


(115) 


and roots 


Xi j2 = 0.625142, 2.7^278 


which can be normalized to 


00i = 22.5051 

CU2 = 99.8205 


Comparing these results with those for the case where simple differences 
were used, one finds that the h 4 approximation gives somewhat better accuracy 
for the first mode and approximately the same results for the higher mode. 

This improvement due to using an h 4 appro ximatic n seems to be rather slight. 
However, if one looks at the accuracy of a given mode as the number of cells 
is increased, the h 4 approximation converges toward the correct result more 
rapidly, as it must, of course, in the limit, fee Table II and Fig. 15 for 
a summary of the results . 

5. Free Beam . — Previously we have seen that the representation of a free 
end is obtained from that for a built-in end by interchanging the roles of y 
and M and also those of 9 and V. Furthermore, ve found for a continuous beam 
that the nonzero mode frequencies are identical for the cases where both ends 
are built-in or both ends are free. This is also true when finite difference 
methods are used. To illustrate this, let us calculate the mode frequencies 
(even symmetry only) of a 6-cell beam that is free at both ends. 

The equations of motion are 
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(lib) 


(119) 


Performing the matrix multiplications for the case where h = 1, we obtain 
from Eq. (8l) that 


(0.206l9-A 2 ) 

-2.16868 

1.23959 


-0.00031 

(5.46771-A 2 ) 

-3-64483 


-0.20588 

-3.29903 

(2.40524-X 2 ) 


0 


( 120 ) 


giving the characteristic equation 


A 2 (A 4 - 8.07914a 2 + 3-00462) = 0 (121) 

Comparing Eqs . (115) and (121), we see that the nonzero roots are identical 
even though the K matrices are quite different. 

All computations in this report applicable to either a built-in or free 
beam were actually set up considering the built-in case, using equations 
similar to (110) to (113). 

4. Cantilever Beam > — One can use the methods developed for obtaining 
h 4 accuracy in the representation of a built-in end or a free end to write the 
equations of motion for a cantilever beam. The minimum number of cells re- 
quired for h 4 accuracy is five. Other methods to be considered later will be 
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applicable with fewer cells. 

As an example of a typical calculation, consider a six-cell cantilever 
beam. The equations of motion are 
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The large numbers in Eqs . (122) and (124) result directly from using the 
one-sided approximations of order h 4 as given in part A of Table I. 


Note the similarities in the matrices of Eqs . (122) and (124). Reading 
from left to right, beginning with the first row of (122), one sees the neg- 
ative of the numbers obtained in reading from right to left beginning with 
the last row of (124). A similar relationship exists between Eqs. (123 ) and 
(125), except for sign. This property is common to all finite difference 
analyses of the cantilever beam and results from the complementary character 
of y and M and also of 9 and V . 


Setting h = 1, one obtains the characteristic equations for the system: 


A 10 - 29.01745A 8 + 221. 62314 A 6 - 462.526o8a 4 + 148.71888a 2 - 1-37737 


0 

( 126 ) 


giving the normalized frequencies 


COi = 3-51687 

CCfe = 22.1688 

C03 = 59-3364 

co 4 = 99.8114 

0*5 = 153-687 


Note that all except the lowest mode frequency are roughly equal to cor- 
responding frequencies for the 6-cell built-in case. An even closer agree- 
ment occurs in the percentage error for corresponding cases. In fact, for 
large N and comparing the higher mode's, one finds essentially identical re- 
sults for the built-in and cantilever beams. (See Tables II and III.) 
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c. ORTHOGONALITY 

We have seen that the use of approximations of order h 2 leads to orthog- 
onal modes. Also, one can obtain a lumped physical system composed of linear , 
bilateral elements that is represented exactly by the equations of motion. 

This is no longer true for the h 4 approximation using one-sided differences 
at the boundaries. If one writes the equations of motion, the m matrix is 
symmetric but the k matrix is not. [See Eq. (120), for example.] Thus the 
system is nonphysical, i.e., a passive system of Linear bilateral elements 
cannot be obtained for which the equations are an exact mathematical repre- 
sentation. Also, the modes are not orthogonal, implying that one cannot ob- 
tain a coordinate transformation that diagonalizes the m and k matrices 
simultaneously. 

To illustrate this point, consider again the 6 -cell built-in beam for 
the case of even symmetry. For each mode, one cai calculate an amplitude 
ratio that indicates the relative magnitude of th ? motions of y 2 and y 3 . 


A£ _ 5.22099-A 2 

As 2.77553 


(127) 


The results for this case are given by the matrix 


[A] 


1.00000 1.00003 

1.74028 -0.88897 


(128) 


where each column corresponds to a natural mode aid the amplitude of the y 2 
motion is arbitrarily set equal to unity. 

The original mass matrix for this system is 
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The deviations from a unit matrix in this case are brought about by the elimina- 
tion of y ± using the constraint equation 


yi = | y 2 - y 3 


(150) 


The stiffness matrix for the original system is 


[k] = [m] [K] = 


5.51698 -2.93800 
-4.34715 2.68740 


( 131 ) 


Consider now the coordinate transformation 


(y) = [A](y] 


(132) 


The mass matrix in terms of the barred coordinate system is 


[m] = [A] [m] [A] 


4.05185 -0.50770 

-0.50770 I.65671 


and the stiffness matrix is 


[E] 


[A] T [k][A] 


1.58347 - 3*90342 

- 0.19842 14.27502 


( 133 ) 


( 134 ) 


Thus we see that neither the mass nor the stiffness matrix is diagonalized 
by the transformation to generalized modal coordinates, and therefore the modes 
are not orthogonal. Furthermore, the k matrix is not symmetric. 

We have defined orthogonality of modes in terms of a simultaneous diagonaliza- 
tion of the mass and stiffness matrices by a coordinate transformation. It should 
be noted, however, that a nonorthogonality of modes in this sense does not nec- 
essarily imply dependence or coupling of modes. In the above case, for example, 
the matrix [!TT]" 1 [Ic] turns out to be diagonal. 
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VI. THE USE OF SYMMETRY IN REPRESENTING BOUNDARY CONDITIONS 


The use of an h 4 approximation in a straighl forward manner has been show r n 
to require the use of one-sided differences near the ends, whereas central 
differences can be used in middle portion. Another approach to the problem 
of representing boundary conditions is to assume a virtual or image beam ex- 
tending beyond the actual beam and to use values of y, 9, M, and V in this 
imaginary beam, as necessary, in the standard central difference equations. 
This approach will now be investigated. 


A. SYMMETRY ASSUMPTIONS 

In establishing the basic approach to the problem, we will again use the 
central difference approximation 


ay 

dX' 




1 

24h 


y n -i - 2Ty n + 2 7y n +i - yi.+ 2 ) + 


2 ) 


3h 4 d 5 y 

640 6 x 5 


|n+|- 


(99) 


using staggered stations such that y and M are calculated at integer stations 
whereas 9 and V are calculated at half -integer slations. The beam ends occur 
at half-integer stations. 


Let us now consider appropriate symmetry assumptions for a cantilever beam, 
chosen to illustrate free and built-in end condil ions . (See Fig. 10.) First 
we note that the assumption of even symmetry in y about the built-in end will 
give zero slope at that end, using Eq. (99)* (The assumption of odd symmetry 
here would not give the proper result.) Successive differentiation of y with 
respect to x gives alternating odd and even fund ions, as shown by the dashed 
lines where they deviate from the solid lines. 'hus the assumed beam shows 
a discontinuity in V (and the slope of M). However, one can see from the basic 
approximation given in Eq. (99) that the presence of discontinuities will 
amount to an increase in the magnitude of the higher derivatives at this point 
and thereby increase the error in the approximation. 

To avoid these errors due to discontinuities , we will arbitrarily assume 
continuous curves at the boundaries even though these assumptions are not 
physically consistent. We are interested in obtaining the best possible 
estimates for the derivatives of y within the actual beam, and therefore we 
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make the most promising assumptions concerning the virtual beam. 

The symmetry assumptions that we will use can be summarized as follows: 

(1) If the derivative of the function approaches zero as the end of the 

actual beam is approached, assume even symmetry about the end. 

(2) If the derivative of the function does not approach zero as the end 

is approached, assume odd symmetry about the end, possibly with an offset to 
avoid discontinuities. 

The procedure by which these assumptions are incorporated into the equa- 
tions of motion is quite straightforward except for the extrapolation of M 
at a built-in end or the similar extrapolation of y at a free end. Taking 
the case of the built-in end, a functional form must be assumed for M near 
the boundary. Referring to the analytical solution near a built-in end, as 
given by Eq. (22) or (49), we find that the magnitude of M varies according 
to a linear plus a fourth-order term in x plus other higher-order terms. So 
a logical choice of the functional form of M near the boundary might be 


M = a Q + a x x ± a 4 x 4 


(155) 


where the last term takes the plus or minus sign depending on whether x is 
positive or negative. This assumption results in the equations 


and 


Mo = ^ [503Mi - 310M2 + 39M 3 ] 


(156; 


M 


252 


[735Mi - 572M 2 + 39M 3 ] 


(137) 


where the built-in end occurs at station 1 .. 

Another possibility is to assume that M is an offset odd function having 
a linear plus cubic variation with x. Thus we might assume that 


M 


b Q + b 3 .x + b 3 x 3 


(138) 



ho 


±n±Q assumption avoids having discontinuous derivatives at x = 0 and gives 
numerically simpler equations. In this case 


l[7Mi - 5 Ms + M 3 ] 

(139) 

1[10Mi - 8M2 + M 3 ] 

(140) 


Computations were run for many cases, using either Eq- (155) or (138)* and 
the results were compared. The differences were fairly small and, if anything, 
favored the linear plus cubic assumption. All results given in this report, 
assuming symmetry at the boundaries, are based upoi a linear plus cubic varia- 
tion of M near a built-in end, i.e., upon Eqs . (13 2) and (140). 

The assumptions concerning the form of y near a built-in end are the 
same as we used previously for the h 4 case. Again we can eliminate y x by 
using the equation 


yi 



(130) 


Of course a similar relation exists for the value of M at the station 
nearest a free end. 


B. EQUATIONS FOR A BUILT-IN BEAM 

1. General Case . — We have seen how we can us b symmetry to extrapolate 
the variables of the problem beyond the actual limits of the beam. Using these 
extrapolated values and the basic central different approximation of Eq. (99 ) > 
we can write general equations for a built-in beam similar to those given by 
h 4 approximation in Eqs. (106) to (109)- They are 
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A comparison of the above equations with those for the h 2 * 4 approximation 
as given by Eqs . (106) to (109) shows considerable similarity. It should be 
noted, however, that the use of symmetry at the boundary never requires more 
than four inputs per amplifier in the analog computer circuit, whereas the h 4 
approximation requires up to five inputs per amplifier. 

2. The 6-Cell Case . — As an example of a mode calculation for a specific 

case, consider again the modes with even symmetry about the center for a 6- 
cell beam built-in at both ends. Th£ equations of motion (using symmetry 

assumptions at the boundaries) are 
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Setting h = 1 and performing the matrix mult [plications , one obtains 


(8. 04779- A 2 ) -4. J9000 

-4.59275 (2. 91454 -A 2 ) 


0 


(149) 


giving the characteristic equation 


A 4 - IO.96255A 2 + 4. 17146 


0 


(150) 
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The roots are 


Ai^ 2 = 0.62828, 3.25078 


and the normalized roots are 


COx = 22.6182 

C£s = 117.028 


Recalling that the corresponding normalized frequencies for the con- 
tinuous beam are 


Pi = 22.3733 

fe = 120.903 


we see that the above results give a somewhat larger error for the first mode 
than does the h 4 approximation, but the error in the higher mode is much smaller. 

An alternate method can be used, particularly with the use of a digital 
computer, to set up the equations of motion and to solve them. This method, 
which indicates more clearly the symmetry assumptions in its formulation, is 
given in Appendix A. 

3. Results . — The results of calculations for built-in and cantilever beams, 
using symmetry assumptions at the boundaries, are summarized in Tables II and 
III and in Figs. 17 and l8. 

Comparing these results with those for the h 4 approximation, we note that 
for a practical range of N (say 6 to lb) the use of symmetry gives lower ac- 
curacy on the low frequency modes but better accuracy on the higher modes. For 
example, consider the case of an 8-cell cantilever beam having a total of 
seven modes. Using symmetry at the boundaries, the first four modes have a 
frequency error well within 1 $ and even the sixth mode has an error of only 10$, 
approximately. By contrast, the use of h 4 approximations at the boundaries 
gives better accuracy for the first three modes but the fourth mode has an 
error of nearly 5$ and the sixth mode has an error of approximately 20$. 
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C. ORTHOGONALITY 

It turns out when boundary conditions are imposed using symmetry that 
the resulting modes of vibration are nearly, but not exactly, orthogonal. 
This means that this approach also results in a nonphysical system in the 
sense that a lumped model cannot be built of linear bilateral elements. 

To compare the orthogonality properties with those for the h 4 approxi- 
mation, let us consider again the example of the 5 -cell, built-in beam. 

The mass matrix is again 


(129) 


The stiffness matrix is 


(151) 


[k] = [m][K] 


8.48426 -4.63270 

-4.47131 2.958231 


85 

-2 

81 

225 

-2 
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225 

625 


The amplitude ratios are calculated from 


As m 8.04779-A 2 

Ai 4.39000 


( 152 ) 


yielding the modal matrix 


[A] 


1.00000 1.00000 
I.7A329 -0.573! >9 


(153) 


In the modal coordinate system the mass matrix is 
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(154) 


(155) 


Thus neither matrix is diagonalized, by the transformation , but the result 
is considerably closer than was the case for the h 4 approximation, as can be 
seen from Eqs . (135) an( l (134)* 
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VII. THE PASSIVE CIRCUIT NETHOD 


We have seen in the case of the simple approximation of order h 2 that an 
electrical circuit exists which the difference equations describe exactly. 
However, when one goes to higher-order differences and uses the methods we 
have described for establishing the boundary conditions, one finds that the 
resulting equations are nonphysical, i.e., no passive linear circuit analogy 
exists . 

It can be shown that the nonphysical character of the higher-order methods 
discussed thus far is due to the boundary-condition representation rather than 
the basic difference approximation. In this section a passive circuit analogy 
will be obtained for the basic higher-order approximations of Eqs . (100) to 
(105). By imposing boundary conditions on this circuit and writing the cor- 
responding equations, one can obtain a set of mutually orthogonal modes. 


A. PASSIVE CIRCUIT ANALOGY 

We have been using the following finite difference approximations to the 
beam equations. 

®n+£ = 2^h (V 1 " 1 " + +1 ’ y n+^) 

M = — (0 o-279 , + 270 ^ - 0 ) 

n 2kh 1+ Z n+ V 

V „ + ! ' ^ (Vl ' 27M n * 2 T">«. - >W) 

Setting h = 1, and representing generalized velocities by voltages and gen- 
eralized forces by currents, one can show that the passive circuit of Fig. 11 
is described by these equations. Note that three transformers, one inductor, 
and one capacitor are required per cell. 

The rather large transformer requirement of Fig. 11 compared with the h 2 
approximation of Fig- 3(b) makes it impractical 1o use this passive circuit 
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directly. Instead we can write the circuit equations and solve them using a 
differential analyzer or a digital computer. The equipment requirements in 
this case are almost the same as for the other higher -order methods hut are 
actually slightly less due to the details of boundary representation. 


B. THE BUILT-IN BEAM 

1. Symmetry Assumptions . —To obtain a proper representation of a built- 
in end, it is convenient to make certain symmetry assumptions consistent 
with the given boundary conditions . We will assume that the deflection curve 
has even symmetry about the end and zero slope at the boundary. Thus y and 
M are even functions of x, while 9 and V are odd functions of x, as shown 
by the dashed lines of Fig. 10. 

The circuit that has been used to represent a built-in end is shown in 
Fig. 12. The beam end is at the center of the diagram (station Jr) with the 
virtual beam on the left and the actual beam on the right. The auto trans- 
former at the bottom of the figure represents the constraint relating y x 
and y^. 


2. Constraint Equation.— It can be shown for a uniform beam that 


y 


n + i 



+ 9y + 9y. 


n+i 



+ 


512 bx e 


n+ 2 


(156) 


So if we assume that y has even symmetry about the built-in end at station i 
and also set yi equal to zero, we obtain the approximate equality 
2 


yi 



(157) 


This is used as the constraint equation in the passive circuit approach. 

It can be seen from Fig. 12 that the autotrans former representing the con- 
straint equation causes a certain relationship to exist between voltages and 
also between currents. In mechanical terms, a constraint equation relating 
coordinates also implies forces of constraint to be exerted on the system in 
such a way that no work is done. These forces of constraint have been ignored 
in the previous two approaches to the problem and this resulted in a lack of 
symmetry in the stiffness matrices. 
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To calculate the effect of these forces of constraint, consider a beam 
described by the matrix equation 


[m](y} = [D ] (V} 


( 158 ) 


Nov suppose a constraint is imposed as described by 


yi = (c){y} 


(159) 


where the summation is over all y except y x . (Tie parentheses denote a row 
matrix. ) Let us include the constraint forces in the analysis and also separate 
the y x equation from the others . Then we obtain 


miyi = 

(Dx)(V} + Fx 

( 160 ) 

M(y) = 

[D](V) + (F) 

(l6l) 


where the F are constraint forces acting on the team. Note that Eq. (l6l) does 
not contain y x and the D matrix does not contain a D x row. 

Now multiply (l6o) Ly y x and premultiply (l6l) hy (y} T . Adding, we obtain 
yi m iyi + (y} T [m](y) = yx(Dx)(V} + (y} T [D]{V) + Fxyx + (y) T (F) (162) 


Since the constraint does no work, we can write 


Fiyi + (y) T {F) = 0 (163) 

and the last two terms of (l 62 ) drop out. Finally, substituting for from 
Eq- (159 ) > dividing out {y}h and rearranging, we obtain 
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m 1 (c){y}(C) T + [m]{y] = (D x ) {V} (C ) T + [D ]{V ) (164) 


This gives the equations of motion with y x omitted and constraint forces in- 
cluded. For example, in the present case with a simple constraint given by 
Eq. ( 157 ), the effective mass at y 2 is 82/81 times its value before the con- 
straint was applied. More complicated constraints will result in inertial 
coupling, but the effective inertia matrix will always be symmetric. 

3. General Case . —The complete equations of motion for a built-in beam 
can be written now using the passive circuit method. They are 
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(168) 




The amount of analog computer equipment needed for the passive circuit 
method is slightly less than for either of the other higher-order methods. 
There are two reasons for this. First, the value of V at a built-in end (or 
of 9 at a free end) need not be calculated. Second, the signs are optimum in 
that no extra amplifiers are required for inversion at any point. 


The equations of motion for stations near 1 free end can be obtained most 
easily from the above equations by exchanging the roles of y and M and also 
0 and V. Of course they could also be obtained directly from the passive cir- 
cuit, remembering the constraint equation relating the values of M at the last 
two stations. 

4. 6-Cell Case . — Again let us consider the even symmetry modes of 

a 6-cell beam built-in at both ends. Using the passive circuit approach, the 
equations of motion are 
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Setting h = 1 , we can solve for the eigenvalues from 


(172) 


(8.47527-A 2 ) -4.54772 

-4.60386 (2-99140-A 2 ) 


0 


(173) 


or 

4 


A 4 - 11. 4666 7A 2 + 4.41583 = 0 


(174) 


The normalized roots are 


<»i = 22.7395 

CDs = 119-765 


In this case the first mode frequency is 1.6$> high and the second mode fre- 
quency is 0.94$> low. 


C . RESULTS 

The results of mode frequency calculations for a uniform beam using the 
passive circuit method are summarized in Tables II and III and in Figs. 19 and 
20 . 


Comparing these results with those obtained previously, using other methods, 
we note that the passive circuit method gives very good over- all accuracy for 
cases where the number of cells is 9 or more. On the other hand, the symmetry 
method appears to give generally better results for smaller values of N. Of 
course, the h 4 method is still best for accuracies of the order of 0.25$> or 
better. 
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One way of comparing the various methods is to obtain an approximate number 
of cells required per wavelength to obtain a given accuracy. For 1 $ error, 
the ordinary h 2 approximation requires 12 cells per wavelength. In contrast, 
the h 4 approximation requires 6-l/2 cells per wavelength, the symmetry meth- 
od requires 5 cells per wavelength, and the passive circuit method requires 
4-l/2 cells per wavelength. 


D. ORTHOGONALITY 

In contrast with the higher-order methods considered previously, the pas- 
sive circuit method produces orthogonal modes. Che orthogonality property is 
a result of the symmetry of the original m and k matrices. 

We can illustrate this point by referring again to the 6 -cell case. From 
Eq. (173) we see that the amplitude ratios for the two modes can be calculated 
using 


Ag. = 8.47527-A 2 

Ai 4.54772 


The modal matrix is 


[A] 


1.00000 1.00000 
1.77590 -0 . 5 r 005 


Ihe original mass matrix is 


[m] 


82 

81 


0 1 


The corresponding stiffness matrix is 


8.57990 -4.30386 

-4.60386 2.99140 


(175) 


(176) 


(177) 


[k] 


(178) 


4=- owe a 



P On VO -4- 


53 


After the transformation 


(y} = [A]{y] 


the resulting mass and stiffness matrices are 


[m] = [ A ] T [m ] [A] 


[E] = [A] T [k][A] 


4.l66l6 0 

0 1.33730 

1.66225 0 

0 14.80080 


(132) 


(179) 


(180) 


The diagonal nature of these matrices implies orthogonality of modes, and also 
that the motion of each of the modal coordinates is independent of the others. 

It should he noted in passing that the matrix 


[A 2 ] = [A]” 1 [K] [A] (181) 

is diagonal for all three methods and the numbers along the main diagonal are 
the squares of the mode frequencies This does not imply orthogonality , 

however, as we have defined it* 
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VIII. SOURCES OF ERROR 


In this section we will consider the principal sources of error in the 
methods that have been suggested for the solution of beam vibration problems. 
These errors are (l) finite difference errors due to the basic difference approx- 
imation; (2) finite difference errors in the representation of boundary con- 
ditions, and (5) component errors. The error curves that have been presented 
so far represent a combination of error sources (l) and (2) above for various 
specific cases. Now we will consider them with somewhat more generality. 


A. FINITE DIFFERENCE ERRORS 

1 . The h^ Approximation . We have seen previously for a pinned beam that, 
when proper symmetry assumptions are made, there is no additional error due to 
the method of representing boundary conditions and so the entire computed 
error is due to the basic finite difference approximation. So consider now 
a uniform beam pinned at each end at an integer station. 

Recall again the basic h 2 approximation 


dx 


n+4 


n+i 


h 


h^ 3fy 

24 3x3 


Repeated application of Eq. (60) yields 


dfy 

dx 2 


n 


'n-i 


2y + y 
J n J n+i 


h 2 3 4 y 
12 dx 4 


Similarly, one can obtain 


* 

(60) 


(65) 



y n - 2 ~ -1 + 6 y n - 4 ^ n+1 + y n+g h g a 6 y 

h 4 6 dx 6 


(182) 


which is the result of applying Eq. (60) four times in sequence. 
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We will again consider the equation for the continuous team in simplified 

form. 


+ afy = 0 
at 2 


As we have seen previously, the solution for the continuous pinned heam is of 
the form 


y = Z Yj_ sin s/ftT x (G^ cos £3j_t + sin f3j_t) (l8j) ) 

i 

where 

Pi = i 2 it 2 (i = 1, 2, 3, •••) (5'J>) 

For the case of the ith mode, we can substitute Eqs • (l82) and (l8d) into (6), 

obtaining 

7 (y n -2 - 4y n-i + 6y n - 4y n+i + y n +2 ) + P i y n ' ^ y n = 0 

or 

^ ( y n- 2 - 4y n-! + 6y n ’ 4y n +1 + y n +2 ) = A y n ' g- P i) < l84) 

Now the left-hand side of Eq. (l84) is of the same form as the approximation 
to av used in the finite difference solution. Therefore, the corresponding 
ax 4 n 

finite difference equation is 

o|y n (183) 
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where 03j_ is the finite difference eigenvalue corresponding to 3-j_ . So the right 
sides of Eqs. (l84) and (185) can he equated if sinusoidal mode shapes are 
assumed in each case. Finally, noting that 


h = 1 

N 


( 186 ) 


we obtain 


03 ? 

1 


(i*r 


1 - uift * + 


] 


or 


0^ = 


t - my * 


(187) 


where i is the mode number. 

4 

2. The h Approximation . — A similar approach can be used to compute the 
frequency error for a pinned beam using higher-order differences. The basic 
approximation is 


8^ 

8x 


n+i 


24h 


27y n + 27y n+1 - / n+2 


2 ) + 

' 64o 8 x 5 


(99) 


n+± 


Using successive repetitions of Eq. (99) * we obtain 


3; 

Sx 2 




( 188 ) 


and 


¥4 

aT 1 


(24h)“ 
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Proceeding in a manner similar to that for the h 2 approximation, we ob- 
tain 


or 


«? 


ft 2 3 h4 r 4 
Pi " *~ 7 ~ Pi 
1 160 



(190) 


Thus we see that the basic finite difference approximation tends to give 
frequencies that are too low. The first error term in Eq. (190) gives an 
accurate estimate of the error only for small i/N ratios and more terms are 
needed for most practical cases. The actual error curves are plotted in Fig. 
21 for the h 4 approximation. 

5. Maximum Errors . — From the tabulated results of the calculations of 
mode frequencies, it is seen that, as the i/N ratio increases, the error also 
increases. In the limit as i/N approaches unity, the error is independent 
of the boundary conditions and depends only on the basic approximation. This 
limiting value of the error could be obtained by evaluating the complete 
series indicated in Eq. (187) or (190) for i/N =1. A simpler way to obtain 
this result is to note that as i/N approaches unity the motion at adjacent 
stations is equal in amplitude but opposite in sign. Knowing this, one can 
calculate the frequency directly. 

For example, using the h approximation and assuming h = 1, we find that 


y n + l6 y n = 0 


(191) 


yielding the frequency 


A 


4 
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The corresponding frequency for the continuous beam is the same as for a pinned 
beam of unit length vibrating in its fundamental mode, namely, n 2 • So the 
maximum error for the h 2 approximation is 


e 


max 



59 - 472 ? 


Performing a similar calculation for the h 4 case, we obtain 


( 192 ) 


y + 

J n 


(i ) 4 


= o 


or 


( 195 ) 


A 


49 

9 


giving a maximum error 


e 


max 


9« 2 -49 

9it 2 


44.836* 


( 194 ) 


B. BOUNDARY-CONDITION ERRORS 

We have seen that an analysis of uniform pirned beams with various numbers 
of cells gives directly the finite difference error due to the basic central 
difference approximation. Furthermore, we note taat the eigenvalues of the 
pinned case alternate with the eigenvalues of eitner the cantilever or built- 
in case. Consequently, if the method of applying boundary conditions intro- 
duced no additional error, one would expect the error curves for the cantilever 
or built-in beam to run between and roughly parallel to the error curves for 
the pinned beam. Instead, one finds that all the boundary-condition methods 
considered here result in raising the natural frequencies. This causes a gen- 
eral improvement in accuracy for i/N values greater than about 0.3 for the 
symmetry and passive circuit approximations. For the h 4 method the accuracy 
is improved for all i/N. 
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For the lower inodes and large N, the error for the three built-in or 
cantilever cases mentioned above exceeds the error for the corresponding pinned 
case because the boundary-condition errors are of lower order in h than the 
basic central difference errors , and ultimately must predominate. For example, 
in the equations for a built-in beam using the simple "h 2 " approximation, we 
find the expression 


e a = Zs. (m) 

2 h 

[See Eq. (82).] The error in this case is actually of order h, as can be ascer- 
tained by the Taylor expansion method. 

Similarly, the passive circuit and symmetry methods of approximating a built- 
in end use the expression 



which has a first error term of order h. On the other hand, each equation used 
in the h 4 approximation is truly of order h 4 , and therefore the boundary- 
condition errors do not tend to dominate for large N. 

A study of the error curves presented in Figs. l6 and 21 liSads one to the 
conclusion that the first error term alone does not give a reliable estimate 
of the actual errors encountered. The i/w ratio must be the order of 0.1 or 
less before the first error term clearly predominates and by this time roundoff 
or component errors mask the desired truncation error. 

To illustrate this point, it is of interest to calculate the rate at 
which the actual truncation errors decrease with increasing N as N approaches 
lo, i.e., at the right-hand edge of the error diagram. Rough calculations 
indicate that the h 2 and symmetry methods have an error that is decreasing 
as h 1 * 9 , -while the error for the passive circuit method is decreasing as h 3 '° 
and that for the h 4 method is decreasing as h 4#7 for increasing N. 

The conclusion from the preceding anal ysis is that the method of repre- 
senting boundary conditions should be chosen on the basis of the require- 
ments of the problem and actual error curves, rather than assuming that an 
h 4 approximation is necessarily more accurate than an h 3 or lower approxima- 
tion over the practical range of N. Our results show T that so-called higher- 



60 


order methods are generally superior to simpler approximation methods. How- 
ever, the choice between the various methods of representing boundary condi- 
tions should be made on the basis of the problem at hand. 


C. COMPONENT ERRORS 

In this section we will consider the sensitivity of the calculated roots 
or frequencies to small shifts or inaccuracies of the component values. This 
will be accomplished by perturbing one of the matrix elements and calculating 
the resulting frequency shift. 

Suppose we write the characteristic equation in the form 


s(n, bj_ j ) - o 


(197) 


where the eigenvalue, r\ = A 2 , is the square of the frequency (rad/sec) of one 

of the modes and b .» 4 is the matrix element that is to be perturbed. For 
r J 

small deviations about the reference values t] 0 , b^j 0 ve can ex P anc ^ i n the 
series 


5s 


SU, Dij) = S(T, 0 , t iJo ) + ^ 


At) + — — 

o 


Ab ij + 


where 


n = n 0 + atj 


b ij 


b ijo + Ab ij 


(198) 


From Eqs. (197) and (198) we obtain an expression for the shift in the root. 


At] 


5s 

1 

Sb ij 

Li 

as 


8ti 

o 


(199) 


To illustrate this method, let us make an error-sensitivity computation 
for a 6-cell beam that is built-in at both ends u£ ing symmetry assumptions at 
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the ends and assuming even symmetry about the center. We have seen that the 
equations of motion can be written in the form 

lyj + [K](y} = 0 (79) 

where 

[K] = - [B 4 ][B 3 ][B 2 ][Bi] (80) 

For the case we are considering, the characteristic equation is of the form 


(Kh-tj) 

K 21 


Kia 


(Kh-tj) 


0 


or 


S = rj 2 - (Kii + ^ 2 . 2 .)^ + (K11K22 ~ K12K21) - 0 


( 200 ) 


Let us suppose that an element in the Bi matrix is to be varied. We can write 


[K] = - [T] [Bi] (201) 

where 

[T] = [B 4 ][B 3 ][B 2 ] (202) 


Consider now the case at hand where 



6 2 


[Bi] 


4775 9 

1 

5400h _ 6q25 ^ 342 


and evaluate the sensitivity to a perturbation of b 21 . Evaluation of At) ac- 
cording to Eq. (199 ) > assuming h = 1, gives 

At)! = 0.157010 Ab 21 

At |2 = 4.22168 Ab 21 


for the first and second modes, respectively- Ir this case the actual roots 
are 


T)1 = 0.394741 

t ) 2 = IO .5676 


So the fractional changes in the roots are 


Mi = _ 0.443791 Msi. 

T)i h 2 i 

^ = - 0.445731 Msa 

This implies that a 1$ increase in the magnitude of h 21 will cause a 0.222$ 
decrease in the first mode frequency and a 0 . 223 $ decrease in the second mode 
frequency. 

o c< 

It has been pointed out by Clymer, Fishery and others that, in general, 
the repeated use of finite difference approximations to the first spatial 
derivative will result in less component error sensitivity than using approx- 
imations to the second or fourth derivatives. Tt check this point, let us 
calculate the error sensitivity for the same problem using approximations to 
the second derivative. Splitting the matrix multiplication differently, we 


p awo ^ 
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will let 


[K] = - [Tj [B J 

where 

[T] = [B 4 ][B 3 ] 

[B] = [Be J [Bo. ] 


In this ease, for h = 1, we find that 


[B] 


r 

58.1271 

-130.2085 

72.0812 


-2.5037 

70.3125 

-67.8088 


x 10 -s 


and 


[TJ 


-1981 4115 -2134 

152 -2182 2030 


Again we will calculate the error sensitivity for b 21 . Performing the 
computations as indicated in the equation, we obtain 


II 

H 

P" 

<3 

- 77.6587 Ab; 

At]2 = 

-4037.34 Ab 21 


corresponding to 


^ = 0.256163 

T)l b 21 


63 

( 203 ) 

( 204 ) 

( 205 ) 


numerical 
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Alls. = 0.497460 AAsi 

Tl2 1321 


So in this case a 1 $ increase in the magnitude oi b 21 will cause a 0.128$ in- 
crease in hfa* first mode frequency and a 0.249$ increase in the second mode 
frequency. 

Comparing these results with those for the approximation to the first 
derivative, we do not have a very significant difference in sensitivity. In 
fact, neither case would be particularly troublesome from the standpoint of 
component error sensitivity. Further calculations of error sensitivities for 
other matrix elements show comparable results. 

Ihis same case, using approximations to the first derivative, was checked 
for error sensitivity on an actual analog computer with results in accordance 
with the above calculations. Percentage frequency shifts were generally less 
than half the percentage component change. 

As the number of cells N is increased, the iifferences in error sensitivity 
between the first derivative and second derivative approximations should be- 
come more apparent. In any case the sensitivity to component errors will in- 
crease with N and large sensitivities should be noticed first in the second 
derivative approximation. 

As Fisher has pointed out, certain elements in the matrix will show con- 
siderably more error sensitivity than others. In general the large terms are 
more sensitive to error. For our basic higher-crder approximation 


& 

6x 


n+ £ 



- 2 ^n + 


= 7y. 


n+i 



(99) 


we would expect a larger sensitivity to errors jn the coefficients ±27 than 
in the coefficients ±1 or even the coefficient 2 k. The coefficients ±1 act as 
correction or trimming terms and do not need to be held to the same percentage 
accuracy. 

The above discussion of error sensitivity applies to digital computers 
when considering roundoff errors. The magnitude of the roundoff error is, of 
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course, dependent upon the number of digits used in the computations and the 
cell number N. It is also dependent upon the differencing scheme used and 
the details of the digital program. For the methods used with IBM 704 com- 
puter in calculating frequencies for the various cases given in this report, 
roundoff errors of the order of O.Ol 1 ^ were first noticed with i/N ratios of 
approximately 0.1 or less on cantilever beam computations. Because the use 
of symmetry effectively halved the number of degrees of freedom, the calcula- 
tions for the built-in beam did not suffer appreciably from roundoff errors. 


For reasons of roundoff error and the general complexity of computation, 
K = 16 seems to be a rough upper limit on the number of cells to be used in 
analog or digital computations. This detracts somewhat from the utility of 
the h 4 method of representing boundary conditions because it is at large N 


that this method is particularly advantageous from the standpoint of finite 
difference or truncation errors. On the other hand, it accentuates the ad- 
vantage of higher -order difference methods in general as compared to the simpler 
h methods, particularly in the computation of the higher frequency modes. 
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IX. CONCLUSIONS 


The principal conclusions gained from the investigations reported here 
are as follows. 

1. The use of higher-order differences in tie solution of beam- vibration 
problems on the analog computer improves accuracy with essentially the 
same amount of equipment as would be required for simpler difference 
methods. In the digital formulation of the corresponding eigenvalue 
problem, a similar improvement in accuracy is obtained with essentially 
the same computer effort, as measured by matrix sizes and numbers of 
operations. The improved accuracy is obi ained by using more inputs 

per amplifier in the analog circuit and more nonzero elements per matrix 
in the digital computations. 

2. The principal difficulty in the formulation of beam-vibration problems 
using higher-order differences lies in the proper representation of. 
boundary conditions. This difficulty arises from the need for obtain- 
ing input data from more stations for each differencing operation to 
obtain greater accuracy than in the simpLer approximations. As the 
boundary is approached, some of these stations needed for the basic 
higher-order central difference approximation lie beyond the end of the 
beam, and so other approximations or assimptions must be used. Three 
methods are presented for the representation of built-in or free ends, 
namely, (l) the use of one-sided differences of order h 3 4 , ( 2 ) the use 
of symmetry assumptions at the boundary, and (3) the passive circuit 
method . 


Each method has its areas of specie 1 utility, as seen from the re- 
suits of calculations for uniform beams. The h 4 method is best suited 
to relatively small i/N ratios, correspcnding to accuracies of 0.25$ or 
better. The symmetry method is most advantageous in obtaining moderate- 
ly good accuracies (the order of 1 $) for all but the highest modes with 
a minimum of equipment, particularly for a cell number N < 9* The pas- 
sive circuit method has quite good accuracy in general, and m particu- 
lar for the higher modes (large i/N ratios) where N > 9. 

3 . The discussion of errors centered on (l y truncation errors from the 

basic central difference approximation told also from the boundary- 
condition representation and (2) component and roundoff errors. Error 
considerations as well as the ease of representing boundaries were im- 
portant in choosing to approximate the :’irst spatial derivative rather 
than the second or higher derivatives and also in the use of staggered 
stations in the calculation of successive derivatives. 
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The methods proposed in this report are relatively insensitive 
to component or roundoff errors for N < 16. Larger values of N are 
not recommended because of the general computational complexity and 
also the increasing size of component or roundoff errors. 

4. The h 4 and symmetry methods are nonphysical in the sense that no system 
composed of linear bilateral elements can be found that is represented 
exactly by the equations of motion. This is a result of the fact that 
the stiffness matrices are not symmetric in these cases. On the other 
hand, the h 2 and the passive circuit methods result in physically 
realizable systems. The latter two methods have the further advantage 
that the resulting modes are orthogonal, i.e., a transformation to 
modal coordinates simultaneously diagonalizes the mass and stiffness 
matrices. 


University of Michigan, 

Ann Arbor, Mich. , October 1961. 
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APPENDIX A 

ALTERNATE MATRIX FORMULA! ION 


The matrix formulation of the difference equations that we have used thus 
far involves the solution of Eq. (8l) for its eigenvalues, the matrix K teing 
obtained after three matrix multiplications. Any assumptions with regard to 
boundary conditions or constraints are put into the individual matrices before 
the matrix multiplications are performed. An alternate scheme, particularly 
adapted to digital computers, enables one to obtain the same K matrix without 
losing sight of the subsidiary assumptions. 

As an example of the alternate method, consider again the problem of the 
6-cell, built-in beam using symmetry assumptions at the ends. Equations (145) 
to (148) could be obtained from the following ei ?ht equations: 
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These equations explicitly separate the basic difference equations and the 
boundary conditions. 
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APPENDIX B 

TAPERED CANTILEVER BEAM 


The analysis that has Been made in the main bcdy of the report has Been 
concerned- with the lateral vibrations of uniform beams. In this appendix we 
will indicate the sort of results to be expected when higher-order methods are 
applied to nonuniform beams. 


ANALYTICAL SOLUTION 

The case of a beam of unit width and with unij orm mass per unit volume and 
a linearly decreasing depth has been analyzed by S.ddall and Isakson. Theoret- 
ical eigenvalues were obtained for the case where 1 he depth at the free end is 
one-fifth that at the built-in end. (See Fig. 22.) For a beam of unit length, 
the linear density is 


p = 2(1 - 0.8x) (Bl) 

and the bending stiffness is 

EI = -(1 - 0.8x) 3 (B2) 

3 


The cantilever beam is built-in at x = 0. 


Theoretical frequencies (rad/sec) for this case are 

= 2.^7829 

p 2 = 9.08902 

3s = 21.2953 


FINITE DIFFERENCE FORMULATION 

The same tapered beam was studied using the finite difference approach. 

In a manner similar to that for the uniform beam, the natural frequencies were 


4^ OWO O 



4=- On V 
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calculated using various boundary-condition assumptions. 

For a nonuniform beam, new assumptions must be made concerning the method 
of lumping the mass and stiffness. Again we assume uniform cell size (h=l) and 
end the beam at half-integer stations. The total mass of each cell is lumped 
at its center. Similarly, the total effective bending stiffness of each cell 
(i.e., between consecutive half-stations) is lumped at its center. 

Specifically, the mass lumped at the nth station is (see Fig. 22) 


“n 


n+ f 

p(x) dx 


n-A 

2 


(B3) 


Using higher-order differences, we obtain the acceleration at the nth station 
from 


y n 


2lm n h 


(- v n-3 + 27V n _ i - 2 ^ V n+i + V n+3) 
2 


(Bl) 


The lumped bending stiffness is 


(El) n 



n-A 


2 


(B5) 


resulting in the bending moment equation 


Mn 


i^(© n _3 - 27© n _i + 279 n+i - 9 n+ a) 
2kh 2 2 2 2 


(b6) 


The slope and shear force equations remain unaltered from the case of a uniform 
beam. 
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In general, the boundary conditions are imposed in the same manner as for 
a uniform beam. An exception occurs in the passive circuit case where the ef- 
fective mass m 2 near the built-in end includes the reflected value (as seen 
from station 2) of the actual mass m 1# A similar 'eflection of bending stiff- 
ness occurs near a free end. 

Also, it should be noted that the constraint equation 


yi = |y 2 - lry 3 (130) 

is of order h 4 for the case of a nonuniform beam instead of order h 6 as was the 
case for the uniform beam. 


RESULTS 

Error curves for the tapered beam are shown in Figs. 23-26. Comparing 
these results with those for the uniform beam, one finds the various methods 
retain many of their characteristics in spite of a general reduction of accura- 
cy. The h 2 method shows good accuracy for the fundamental mode while the higher- 
order methods (the symmetry and passive circuit melhods, in particular) have 
better accuracy on the higher modes. The h 4 methoc suffers the largest loss 
of accuracy due the nonuniform mass and stiffness c istributions . 

The general lowering of mode frequencies compared to the uniform beam is 
probably due primarily to the mass lumping method which results in a slight 
positive shift in the center of mass location compered to the continuous beam. 
This is because the mass is lumped at the center oi the cell whereas the actual 
center of mass for the cell is located slightly inboard of the center, i.e,, 
toward the built-in end. Similarly, the compliance (inverse stiffness) for a 
cell is lumped at the center and is thus somewhat inboard of its compliance 
centroid. This also results in a lowering of the r atural frequencies. 


o o\\o -d- 
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TABLE I 

ONE-SIDED APPROXIMATIONS OF ORDER h 4 


A. General Case 


= - 7 -[- 93yi + 229y 2 
Bxl^ 24h 
2 


225y 3 + ll.'.y 4 - 22y 5 ] + ^h 4 

640 Bx 5 i x 


By 

Bx 


24 h 


[- 22y x + 17y 2 + 9y 3 - 5y 4 + J 5 ] ~ — ^-h 4 

1920 Bx 5 


Be 

Bx 


rfr[- 229 1 + 179 3 + 9% 

£i 4 n 2 22 


5 ©t + 9 g] 

2 2 


7i_ h 4 B e 


1920 Bx 5 


B. 


For y A 

2 


By | 

Bx i A 

2 


840 h 


[ 3675 y 1 - 1225 y 2 + 44 ly 3 - 


75y 4 ] + ~ ~h 4 

128 Bx 5 


By j 


_L _ ^ 4 S 5 y 


. , = t— — [- I085yi + I0l5y 2 - 63 y 3 + 5y 4 ] + — — h' 

Bx i 3 84 oh 1920 Bx 5 

2 


C. For y^ = 0, 9_i 

2 2 


= 0 


y i = 9 y2 


1 . _L h 4 B 4 y 

25 y3 48 h Bx* 


+ _ 5 _h 5 + J^L h 6 ^ 

_i 160 Bx 5 j x 2880 Bx 6 
2 2 
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TABLE I (Continued) 


D. 


For = 0 

2 


s© 

Sx 


1 


— [- 91© i + 99©3 - 9©s + ©zl 
72L 2 2 2 2 


_L_h- 

1920 


d 5 e 
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TABLE II 

MODE FREQUENCIES OF A BUILT-IN OR FREE UNIFORM BEAM* 






Appro xima :ion Method 


N 

Mode 

Continuous 

h 2 

h 4 

Symmetry 

Passive 

Circuit 

3 

1 

22.373 

- 1 . 466 $ 
22.045 


3.831$ 

23.230 

13.565$ 

25.408 

4 

1 

22.373 

1 . 135 $ 

22.627 


1.850$ 
22 . 787 

4.863$ 

23.461 

4 

2 

61.673 

-17.961$ 

50.596 


- 1.516$ 
60.738 

1.527$ 

62.615 

5 

1 

22.373 

1.256$ 

22.654 

0.796$ 

22.551 

1.345$ 

22 . 674 

2.700$ 

22.977 

5 

2 

61.673 

- 9 - 357 $ 
55.902 

-11.052$ 

54.857 

0.028$ 

61.690 

3.259$ 
63 . 683 

5 

3 

120.90 

-27.833$ 

87.246 

-24.607$ 

91-155 

- 8.313$ 
110.85 

- 6.907$ 
112.55 

6 

1 

22.373 

1 . 059 $ 

22.610 

0.590$ 

22.505 

1.095$ 

22 . 618 

1.636$ 

22.739 

6 

2 

61.673 

- 5.609 $ 

58.214 

- 3 - 755 $ 
59-357 

0.328$ 

61.875 

2.679$ 

63.325 

6 

3 

120.90 

-17.883$ 

99.281 

- 17 . 436 $ 

99.820 

- 3.201$ 
117.03 

- 0.935$ 
119.77 

6 

4 

199.86 

- 34 . 099 $ 

131.71 

-23.131$ 

153.63 

- 14 . 000 $ 

171.88 

-13.274$ 

173.33 

*The 

percentages refer to 

the error in 

frequency as 

compared to the 

continuous 


beam. 


Q (J\\D 
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TABLE II (Continued) 






Approximation 

Method 


N 

Mode 

Continuous 

h 2 

h 4 

Symmetry 

Passive 

Circuit 

7 

1 

22.373 

0.858$ 

22.565 

0.540$ 

22.449 

0.912# 

22.577 

1 . 064 # 
22 . 6 ll 

7 

2 

61.675 

- 3.707$ 
59.387 

- 0.751# 
61.210 

o. 46 o# 

61.957 

2.006# 

62.910 

7 

3 

120 . 90 

-12.333$ 

105.99 

- 9.992# 
108.82 

- 1.472# 
119.12 

0 . 711 # 

121.76 

7 

4 

199.86 

-24.507$ 

151.28 

-19.494$ 
160. 90 

- 6.915# 

186 . 04 

- 5 . 424 $ 
189.02 

7 

5 

298.56 

-38.554$ 

184.05 

-25.617$ 

228.05 

-18. 4 12# 
245.59 

-18.000# 

244.82 

8 

1 

22.373 

0.697$ 

22.529 

0.192$ 

22 . 4 l 6 

0 . 764 $ 

22.544 

0 . 729 $ 

22.556 

8 

2 

61.673 

- 2.627$ 
60.055 

0 . 091 # 

61.729 

0.522# 

61.995 

1 . 493 # 

62.594 

8 

3 

120.90 

- 8.985 $ 
110.04 

- 4 . 806 $ 
115.09 

- 0.695# 
120.06 

1.092# 

122.22 

8 

4 

199.86 

-18.065$ 

165.76 

-l 4 .l 4 o$ 

171.60 

- 5.825# 
192.22 

- 2.086# 
195.69 

8 

5 

298.56 

-29.165$ 

211.49 

- 20 . 117 $ 

258.50 

- 10 . 480 $ 

267.27 

- 9 . 489 $ 
270.23 

8 

6 

416.99 

- 41 . 4 ll$ 
244.51 

- 24 . 953 $ 

313.02 

-21.825# 

525.99 

-21.569$ 

327.05 
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TABLE II ( Continued j 
Approximation Method 


N 

Mode 

Continuous 


h 4 

Symmetry 

Passive 

Circuit 

9 

1 

22.373 

0.572$ 

22.501 

0.10'$ 

22.39' 

0.653$ 

22.519 

0.523$ 

22.490 

9 

2 

61.673 

- 1 - 959 $ 
60.465 

0 . 24 : $ 

61.82: 

0.537$ 

62.004 

1.129$ 

62.369 

9 

3 

120.90 

- 6.832$ 
112.64 

- 2 . 02( $ 
118.45 

- 0.281$ 

120.56 

1.092$ 

122.22 

9 

4 

199.86 

-13.905$ 

172.07 

- 9 . 03 J$ 
l 8 l. 8 l 

- 2.272$ 
195.32 

- 0.620$ 
198.62 

9 

5 

298.56 

-22.773$ 

230.57 

- 16 . 342 $ 

249.77 

- 6.391$ 
279.48 

- 5.078$ 
283.40 

9 

6 

416.99 

-32.907$ 

279.77 

-20.665$ 

330.82 

-13.669$ 

359.99 

-12.996$ 

362.80 

9 

7 

555.17 

-43.707$ 

312.52 

- 26 . 471 $ 

4 o 8.21 

- 24 . 504 $ 

719.13 

-24.337$ 

420.06 


10 

1 

22.373 

0.474$ 

22.479 

0.063$ 

22.387 

0 . 550 $ 

22.496 

0.384$ 

22.459 

10 

2 

61.673 

- 1.518$ 
60.737 

0.217 6 
61.807 

0 . 525 $ 

61.997 

0 . 866$ 
62.207 

10 

3 

120.90 

- 5.360$ 

114.41 

- 0 . 753 $ 
119.99 

- 0.033$ 
120.86 

0.976$ 

122.08 

10 

4 

199.86 

-11.018$ 

177.84 

- 5 . 244 $ 
189.38 

- 1 . 401 $ 
197.06 

0.030$ 

199.92 

10 

5 

298.56 

- 18 . 214 $ 

244.18 

-12.165$ 

262.24 

- 4 . 123 $ 
286.25 

- 2.726$ 
290.42 

10 

6 

416.99 

-26.629$ 

305.95 

-17.506$ 

343.99 

- 8.947$ 

379-68 

- 7 . 959 $ 
383.80 

10 

7 

555.17 

-35.865$ 

356.06 

-21.478?;, 

735.93 

- 16 . 449 $ 

463.85 

- 15 . 972 $ 

466.50 

10 

8 

713.08 

-45.493$ 

388.68 

- 27 . 959 $ 

513.71 

- 26 . 649 $ 

523.05 

-26.534$ 

523.87 


Q On vo 
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TABLE II (Continued) 


N Mode 


11 1 

11 2 

11 3 

11 4 

11 5 

11 6 

11 7 

11 8 




Approximation Method 


Continuous 

h2 

h4 

Symmetry 

Passive 

Circuit 

22.373 

0.402$ 

22.463 

0.035$ 

22.381 

0.478$ 

22.480 

0.295$ 

22.439 

61.673 

- 1.213 $ 

60.925 

0.162$ 

61.773 

0.499$ 

61.981 

0.678$ 
62 . 091 

120.90 

- 4.334$ 
115.66 

- 0.232$ 
120.62 

0.116$ 

121.04 

0.844$ 

121.92 

199.86 

- 8.936$ 
182 . 00 

- 2.867$ 
194.13 

- 0.866$ 
196.13 

0.314$ 

200.49 

298.56 

-14.871$ 

254.16 

- 8.377$ 
273.55 

- 2.767 $ 
290.30 

- 1.440$ 
294.26 

416.99 

-21.919$ 

325.59 

-14 . 185$ 
357.84 

- 6.098$ 
391.56 

- 4.959$ 
396.31 

555.17 

-29.812$ 

389.66 

-18.328$ 

453.42 

-11.366$ 

492.07 

-10.620$ 

496.21 

713.08 

-381248$ 

440.34 

-22.501$ 

552.63 

-18.844$ 

478.71 

-18.499$ 

581.17 

890.73 

-46.918$ 

472.82 

-29.311$ 

629.65 

-28.398$ 

637.78 

-28.317$ 

638.50 


11 


9 
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TABLE II (Continued) 


N Mode 


12 1 

12 2 

12 3 

12 4 

12 5 

12 6 

12 7 

12 8 

12 9 


Appro <imation Method 


Continuous 

h 2 

h 4 

Symmetry 

Passive 

Circuit 

22.373 

0.340$ 

22.449 

0.02?$ 

22.373 

0 . 4 lb$ 

22.466 

0.228$ 

22.424 

61.673 

- 0.992$ 
61.061 

o.n+$ 

61.743 

0 . 467 $ 

61.961 

0 . 538 $ 

62.005 

120.90 

- 3.573$ 

116.58 

- 0.025$ 
120.87 

0 . 199 $ 

121.14 

0.720$ 

121.77 

199.86 

- 7.395$ 

185.08 

- 1.516$ 
196.83 

- 0.520$ 
198.82 

0 . 430 $ 

200.72 

298.56 

- 12.363$ 
261.65 

- 5 . 446 $ 
282.3c 

- 1.906$ 
292.87 

- 0 . 717 $ 
296.42 

416.99 

-18.322$ 

340.59 

-10.856$ 

371.72 

- 4 . 295 $ 

399.08 

- 3 . 132 $ 

403-93 

555.17 

-25.097$ 

415.84 

- 15 . 454 $ 

469.15 

- 8.086$ 
510.28 

- 7 . 171 $ 
515.36 

713.08 

- 32 . 468 $ 

481.56 

- 19 . 124 $ 

576.71 

- 13 . 590 $ 

616.17 

-13.021$ 

620.23 

890.73 

- 40 . 206 $ 

532.60 

-23.639$ 

680.35 

-20. 910$ 
704.48 

-20.653$ 

706.77 

1088.1 

- 48 . 084 $ 

564.93 

-30.5c 8$ 
756 . 1 C 

-29.850$ 

763.32 

-29.786$ 

763.98 


12 


10 
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TABLE II (Continued) 






Approximation Method 


w 

Mode 

Continuous 

h 2 

h 4 

Symmetry 

Passive 

Circuit 

13 

1 

22.373 

0.295# 

22.439 

0.013# 

22.376 

0.362# 

22.454 

0.179# 

22.413 

13 

2 

61.673 

- 0.827# 
61.163 

0.078# 

61.721 

0.431# 

61.939 

0 . 435 # 

61.941 

13 

3 

120.90 

- 2.994# 
117.28 

0 . 04 l# 

120.95 

0 . 248 # 

121.20 

0.6o4# 

121.63 

13 

4 

199.86 

- 6.219# 

187.43 

- 0.791# 

198.28 

- 0.285# 
199.29 

0 . 460 # 

200.78 

13 

5 

298.56 

-10.430# 

267.42 

- 3.^23# 
288.34 

- 1 . 330 # 
294.59 

- 0.301# 
297.66 

13 

6 

416.99 

-15.526# 

352.25 

- 7.897# 
384.06 

- 3 . 103 # 
4 o 4 . 05 

- 1.990# 
408.69 

13 

7 

555.17 

-21.379# 

436.48 

-12.656# 

484.91 

- 5 - 904 # 
522.39 

- 4 . 914 # 

527.89 

13 

8 

713.08 

-27.830# 

514.63 

-16.455# 

595.74 

- 10 . 014 # 

641.67 

- 9.278# 

646.92 

13 

9 

890.73 

-34.710# 

581.56 

-19.988# 

712.69 

-15 . 612# 
751.67 

- 15 . 170 # 

755.61 

13 

10 

1088.1 

- 4 l. 84 l# 

632.83 

- 24 . 740 # 

818.86 

-22 . 700# 
841.13 

-22.498# 

843.26 

13 

11 

1305.3 

- 49 . 054 # 
665 . 02 

- 31 . 564 # 

893.32 

- 31 . 073 # 

899.71 

-31.027# 

900.30 



82 


TABLE II (Continued) 


N Mode 


14 1 

14 2 

14 3 

14 4 

14 5 

14 6 

14 7 

14 8 

14 9 


14 10 

14 11 




Approximation 

Method 


Continuous 

h 2 

h* 

Symmetry 

Passive 

Circuit 

22.373 

0.255$ 

22.430 

0.(09$ 

22.;-75 

0.322$ 

22.445 

0.147$ 

22.4 o 6 

61.673 

- 0.700$ 
6 i .24 i 

0.052$ 
61. "05 

0.397$ 

61.918 

0.355$ 

61.892 

120.90 

- 2.548$ 
117.82 

0.(158$ 

120.97 

0.281$ 

121.24 

0.521$ 

121.53 

199-86 

- 5.309$ 
189.25 

- 0.1 10$ 
199. (>4 

- 0.130$ 
199.60 

0.440$ 

200.76 

298.56 

- 8.919$ 
271.93 

- 2.124$ 
292.^2 

- 0.931$ 
295-79 

- 0.060 $ 
298.38 

416.99 

-13.314$ 

361.47 

- 5-^37$ 
393.90 

- 2.281$ 
407.48 

- 1.266$ 
411.71 

555-17 

-18.407$ 

452.98 

- 9.936$ 
500.01 

- 4.408$ 
530.70 

- 3.412$ 
536.23 

713.08 

-24 . 074$ 
541.41 

-13.955$ 

613-57 

- 7.531$ 
659.38 

- 6.696$ 
665.33 

890.73 

-30.191$ 

621.81 

-17.294$ 

736.09 

-11.844$ 

785.23 

-11.249$ 

790.53 

1088.1 

-36.623$ 
689. 63 

-20.908$ 

860.01 

-17.434$ 

898.40 

-17.085$ 

902.19 

1305.3 

-43.231$ 

741.04 

-25.837$ 

968.06 

-24.263$ 

988.64 

-24.107$ 

990.63 

1542.1 

-49.867$ 

773-10 

-32.482$ 

1041.2 

-32.112$ 

1046.9 

-32.073$ 

1047.5 


14 


12 
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TABLE II (Continued) 






Approximation Method 


N 

Mode 

Continuous 

h a 

h 4 

Symmetry 

Passive 

Circuit 

15 

1 

22.373 

0.223$ 

22.423 

0. 000$ 
22.373 

0.282$ 

22.436 

0.116$ 

22.399 

15 

2 

61. 673 

- 0.602$ 
61.302 

0.034$ 
61. 694 

0.365$ 

61.898 

0.293$ 

61.854 

15 

3 

120.90 

- 2.200$ 
118.24 

0.058$ 

120.97 

0.289$ 

121.25 

0.438$ 

121.43 

15 

4 

199-86 

- 4.578 $ 
190. 71 

- 0.215$ 
199.43 

- 0.020$ 
199.82 

0.425$ 

200.71 

15 

5 

298. 56 

- 7.714$ 
275-53 

- 1.320$ 
294.62 

- 0. 646$ 
296.63 

0.077$ 
298. 79 

15 

6 

416.99 

-11.540$ 

368.87 

- 3.801$ 
401.14 

- 1.698$ 
409.91 

- 0.794$ 
413.68 

15 

7 

555-17 

-16.001$ 

466.34 

- 7.536$ 
513.33 

- 3.347$ 
536.59 

- 2.392$ 
541.89 

15 

8 

713.08 

-21.005$ 

563.30 

-11.526$ 

630.89 

- 5.765$ 
671.97 

- 4.891$ 
678. 20 

15 

9 

890.73 

-26. 455$ 
655.09 

-14.961$ 

757.47 

- 9-128$ 
809.42 

- 8.429$ 

815.65 

15 

10 

1088.1 

-32.249$ 

737.24 

-18.114$ 

891.04 

-13.556$ 

940.62 

-13.069$ 

945.89 

15 

11 

1305.3 

-38.275$ 

805.74 

-21.857$ 

1020.0 

-19.084$ 

1056.2 

-18.800$ 

1059.9 

15 

12 

1542.1 

-44.413$ 

857.23 

-26.859$ 

1127.9 

-25.621$ 

1147.0 

-25.498$ 

1148.9 

15 

13 

1798.7 

-50.564$ 

889.17 

-33.296$ 

1199.8 

-33-007$ 

1205.0 

-32.979$ 

1205.5 
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Approximation Method 


N 

Mode 

Continuous 

h 2 

h 4 

Symmetry 

Passive 

Circuit 

16 

1 

22.373 

0.197$ 

22.417 

0.001$ 

22.375 

0.250$ 

22.429 

0.098$ 

22.395 

16 

2 

61.673 

- 0.524$ 
61.350 

0.025$ 

61.687 

0.334$ 

61.879 

0.245$ 

61.824 

16 

3 

120.90 

- 1.9H$ 
118.59 

0.04 1$ 
120.95 

0. 289$ 
121.25 

0. 380$ 
121.36 

16 

1* 

199.86 

- 3.993$ 
191.88 

- 0.110$ 
199.64 

0.050$ 

199.96 

0.390$ 

200.64 

16 

5 

298. 56 

- 6.736$ 
278.45 

- 0.827$ 
296.OS 

- 0.439$ 
297.25 

0.151$ 

299.01 

16 

6 

416.99 

-10.096$ 

374.89 

- 2.583$ 
406.2C 

- 1.273$ 

411.68 

- 0.484$ 
4l4.97 

16 

7 

555.17 

-14.030$ 

477.28 

- 5.575$ 

524.22 

- 2.576$ 
540.87 

- 1.690$ 
545.79 

16 

8 

713.08 

-18.471$ 

581.37 

- 9.254$ 
647.05 

- 4.483$ 

68l.ll 

- 3.612 $ 
687.32 

16 

9 

890.73 

-23.343$ 

682.81 

-12.761$ 

777.06 

- 7.141$ 
827. 12 

- 6.382$ 
833.88 

16 

10 

1088.1 

-28.562$ 
777. 32 

-15.817$ 

915.96 

-10.671$ 

971.99 

-10.082$ 

978.39 

16 

n 

1305.3 

-34.046$ 
860. 86 

-18.9116$ 

1058.0 

-15.146$ 

1107.6 

-14.748$ 

1112.8 

l6 

12 

15U2.1 

-39.700$ 

929.88 

-22. 8C0$ 
1190. 5 

-20.556$ 

1225.1 

-20.329$ 
1228. 6 

l6 

13 

1798.7 

-45.438$ 

981.40 

-27.8c 9$ 
1298. 5 

-26.825$ 

1316.2 

-26.725$ 

1318.0 

16 

14 

2075.1 

-51.173$ 

1013.2 

-34.018$ 

1369.2 

-33.786$ 

1374.0 

-33.762$ 

1374.5 
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TABLE III 

MODE FREQUENCIES OF A UNIFORM CANTILEVER BEAM 






Approximation Method 


N 

Mode 

Continuous 

h 2 

h 4 

Symmetry 

Passive 

Circuit 

2 

1 

3.5160 

13.766$ 

4.0000 


3.851$ 

3.6514 

4.209$ 

3.6640 

3 

1 

3.5160 

6.030$ 

3.728 


1.812$ 

3.5797 

1 . 257 $ 

3.5602 

3 

2 

22.031* 

- 1.389$ 
21.728 


3.980$ 

22.911 

10 . 674 $ 

24.386 

k 

1 

3.5160 

3.356$ 

3.634 


1.135$ 

3.5559 

0 . 549 $ 

3.5353 

k 

2 

22.034 

1.035$ 

22.262 


1.938$ 

22.461 

5 . 328 $ 

23.208 

h 

3 

61.697 

-17.941$ 

50.628 


- 1.520$ 

60.759 

1 . 483 $ 

62.612 

5 

1 

3.5160 

2 . 142 $ 

3.5913 

0. 060$ 
3.5181 

0.774$ 

3.5^32 

0 . 284 $ 

3.5260 

5 

2 

22.034 

1 . 153 $ 

22.288 

0.880$ 

22.228 

1 . 402 $ 

22.343 

2.927$ 

22.679 

5 

3 

61.697 

- 9 . 333 $ 
55.939 

-11.665$ 

54.500 

0.023$ 
61 . 711 

3 . 232 $ 

63 . 691 

5 

4 

120.90 

- 27 . 841 $ 
87. 242 

-23.956$ 

91.937 

- 8 . 313 $ 
IIO.85 

- 6.898$ 

112.56 
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TABLE III (Continued) 


App roximation Method 


N 

Mode 

Continuous 

h 2 

h 4 

Symmetry 

Passive 

Circuit 

6 

i 

3.5160 

1.1+818$ 

3.5681 

0.028$ 

3 .M 70 

0.555$ 

3.5355 

0.168$ 

3.5219 

6 

2 

22.034 

0.976$ 

22.21+9 

0. 613$ 
22.169 

1.135$ 

22.284 

1.770$ 

22.424 

6 

3 

61 . 697 

- 5.589$ 
58.249 

- 3.827$ 

59.336 

0 . 324 $ 

61.897 

2 . 663 $ 

63 . 3^0 

6 

4 

120.90 

-17.885$ 

99.277 

- 17 . j - 43 $ 

99.811 

- 3 . 201 $ 
117.03 

- 0 . 935 $ 
119.77 

6 

5 

199.86 

-34.099$ 

131.71 

-23.101$ 

153.69 

- 14 . 000 $ 

171.88 

- 13 . 274 $ 

173-33 


7 1 

7 2 


7 3 

7 4 


7 


3.5160 

1.086$ 

3-5542 

0.017$ 

3.1166 

0.4l8$ 

3.5307 

0.102$ 

3.5196 

22.034 

0.794$ 

22.209 

0.545$ 

22..10 

0.944$ 

22.242 

1.148$ 

22.287 

6l. 697 

- 3.691$ 
59.420 

- 0. (72$ 
61.221 

0.459$ 

61.980 

1.997$ 

62.929 

120.90 

-12.333$ 

105.99 

- 9.792$ 
108. 32 

- 1.472$ 
119.12 

0.711$ 

121.76 

199.86 

-24.307$ 

151.28 

-19. *79$ 
160. J3 

- 6.915$ 
186.04 

- 5.424$ 
189.02 

298.56 

-38.354$ 
184. 05 

-23. 520$ 
228. 14 

-18. 412$ 

243.59 

-18.000$ 

244.82 


7 
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TABLE III (Continued) 






Approximation 

l Method 


N 

Mode 

Continuous 

h 2 

h 4 

Symmet ry 

Passive 

Circuit 

8 

1 

3.5160 

0.828$ 

3.5451 

0. 009$ 

3.5163 

0.319$ 

3.5272 

0.063$ 

3.5182 

8 

2 

22.031+ 

0. 644$ 
22.176 

0.191$ 

22.076 

0. 790$ 
22.208 

0.785$ 

22.207 

8 

3 

61.697 

- 2. 6l4$ 
60.084 

0. 086$ 
61.750 

0.520$ 

62.018 

1.488$ 

62.615 

8 

k 

120.90 

- 8.991$ 
110. 03 

- 4.797$ 
115.10 

- 0.695$ 
120. 06 

1.092$ 

122.22 

8 

5 

199.86 

-18.063$ 

163.76 

-14.140$ 

171.60 

- 3.823$ 
192.22 

- 2.086$ 
195.69 

8 

6 

298.56 

-29.163$ 

211.49 

-20.117$ 

238.50 

-10.480$ 

267.27 

- 9.489$ 
270.23 

8 

7 

416.99 

-4l. 4ll$ 
244.31 

-24.931$ 

313.03 

-21.823$ 

325.99 

-21.569$ 

327.05 

9 

1 

3.5160 

0. 648$ 

3.5388 

0.006$ 

3.5162 

0.236$ 

3.5243 

0.043$ 

3.5175 

9 

2 

22.034 

0.531$ 

22.151 

0. 104$ 
22.057 

0. 667$ 
22.181 

0.558$ 

22.157 

9 

3 

61.697 

- 1.948 
60.495 

0.242$ 

61.846 

0. 534$ 
62.028 

1.130$ 

62.394 

9 

4 

120.90 

- 6.832$ 

112.64 

- 2.026$ 
118. 45 

- 0.281$ 
120.56 

1.092$ 

122.22 

9 

5 

199.86 

-13.905$ 

172.07 

- 9.031$ 
181.81 

- 2.272 $ 

195.32 

- 0.620$ 
198. 62 

9 

6 

298.56 

-22.773$ 

230.57 

-16.342$ 
249. 77 

- 6.391$ 

279-48 

- 5.074$ 
283.41 

9 

7 

416.99 

-32.907$ 

279.77 

-20.665$ 

330.82 

-13.669$ 

359-99 

-12.996$ 

362.80 

9 

8- 

555.17 

-43.707$ 

312.52 

-26. 471$ 
4o8.2l 

-24.504$ 

419.13 

-24.338$ 

420.05 
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TABLE III (Continued) 


N Mode 


10 1 

10 2 

10 3 

10 4 

10 5 

io 6 

10 7 

10 8 


Approximation Method 


Continuous 

h 2 

h 4 

Symmetry 

Passive 

Circuit 

3.5160 

0. 509$ 
3.5339 

0 . OC 3 $ 
3.5161 

0 . 148 $ 
3.5212 

0.028$ 

3.5170 

22.031* 

0. 1*1*0$ 
22.131 

0 . 068 $ 

22.049 

0.567$ 

22.159 

0.409$ 

22.124 

61.697 

- 1.509$ 
60. 766 

0 . 217 $ 

61.831 

0 . 525 $ 

62.021 

0 . 864 $ 

62.230 

120.90 

- 5.368$ 
111*. 1*1 

- 0 . 753 $ 
119.95 

- 0.033 $ 
120.86 

0.976$ 

122.08 

199.86 

-11.018$ 

177 . 81 * 

- 5 . 244 $ 

189.38 

- 1 . 401 $ 
197.06 

0.030$ 

199.92 

298.56 

-18. 214 $ 
244.18 

-12.165$ 

262.24 

- 4 . 123 $ 
286.25 

- 2.726$ 
290.42 

1*16.99 

-26.627$ 

305.96 

- 17 . 5 C 5 $ 

343.99 

- 8 . 947 $ 
379.68 

- 7 . 959 $ 
383.80 

555.17 

- 35 . 865 $ 

356.06 

- 21 . 473 $ 

435.93 

- 16 . 447 $ 

463.86 

- 15 . 972 $ 

466.50 

713.08 

-45.491$ 
388. 69 

- 27 . 95 ?$ 
513.71 

- 26 . 649 $ 

523.05 

-26.534$ 

523.87 
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TABLE III (Concluded) 


N Mode 


11 1 
11 2 


11 3 

11 4 

11 5 

11 6 

ll 7 

11 8 


ll 


9 


Approximation Method 


Continuous 

h 2 

h 4 

Symmetry 

Passive 

Circuit 

3.5160 

0.395$ 

3.5299 

0.003$ 

3.5161 

0.054$ 

3.5179 

0.023$ 

3.5168 

22.034 

0.368$ 

22.115 

0.050$ 

22.045 

0. 486$ 

22. l4l 

0.309$ 

22.102 

61.697 

- 1.204$ 
60.954 

0.162$ 

61.797 

0.499$ 

62.005 

0. 676$ 

62.114 

120.90 

- 4.334$ 
115.66 

- 0.232$ 
120.62 

0. 116$ 

121.04 

0.835$ 

121.91 

199.86 

- 8.936$ 
182.00 

- 2.867$ 
194.13 

- 0 . 866$ 
198.13 

0.315$ 

200.49 

298.56 

-14.871$ 

254.16 

- 8.377$ 
273.55 

- 2.767$ 
290.30 

- 1.444$ 
294.25 

416.99 

-21.919$ 

325.59 

-14.185$ 

357.84 

- 6.096$ 
391.57 

- 4.959$ 
396.31 

555.17 

-29. 811$ 
389.67 

-18.328$ 

453.42 

-11.366$ 

492.07 

-10.618$ 

496.22 

713.08 

-38.250$ 

440.33 

-22.501$ 

552.63 

-18.844$ 

578.71 

-18.497$ 

581.18 

890.73 

-46.918$ 

472.82 

-29.311$ 

629.65 

-28.398$ 

637.78 

-28.316$ 

638.51 


11 
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o. Mechanical representation of the simple finite 
difference approximation to a uniform beam. 


^n- 1/2 ®n+i/2 



b. Electrical representation of the simple finite 
difference approximation to a uniform beam. 


Fig. 3. 
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o. Even symmetry about center (symmetric mode). 



i 


b. Odd symmetry about center (antisymmetric mode). 


Fig. 4. 
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a. Analog computer circuit for a four-cell 

2 

cantilever beam using differences of order h 
and second-order equations. 


Fig. 5. 


Analog computer circuit for a f our-cel! cantilever 
beam using differences of order h £ and first- 
order equations. The built-in end occurs at 
station 1/2 and the free end at station 9/2. 
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R 



b. Odd symmetry 


Fig. 6. Analog computer circuit for a six-cell beam built-in at both 
ends. Differences of order h* ore used. Only the left half of the 
beam is represented. 
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Fig. 10. Symmetry assumptions for a cantilever beam. Solit 
lines show assumptions used in computations. Dashed lines 
indicate deviations giving a physically consistent set of 



assumptions. 
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Fig. II. Passive circuit analogy for the basic h 4 approximation. 
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Fig 12. Passive circuit analogy for a built-in end at station 1/2. 
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Fig 13. Built-in, h 2 approximation. 
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F'9 23 Tapered cantilever, h 2 approximation. 
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Fig. 24. Tapered cantilever, h 4 approximation. 
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